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ABSTRACT 


Several  nethods  were  used  to  predict  the  thermal  expansion  coefficients 
of  pure,  single-phase  ceramics.  The  predictions  were  extended  to  cover  1700 
pure  phases.  Based  upon  these  calculations  a  number  of  phases,  predicted 
to  have  low  values  of  thermal  expansion  coefficient,  were  synthesized  and  the 
thermal  expansion  properties  of  these  materials  were  measured.  Several  of 
these  were  fbund  to  be  low-expansion  materials.  The  thermal  expansion 
properties  of  cubic  • OPgO^  are  particularly  interesting.  This  material  expands 
with  increasing  temperature  up  to  about  U00°C.  Above  this  temperature  the 
crystal  contracts,  returning  to  its  room  temperature  length  at  about  1000° C. 
These  properties  are  believed  to  be  unique. 

Variation  of  the  thermal  expansion  anisotropy  by  addition  of  solid 
solution  atoms  to  several  base  crystals  was  attempted.  Addition  of  vanadium 
to  rutile  (TiO^)  resulted  in  a  significant  decrease  in  thermal  expansion 
anisotropy.  Vanadium  additions  to  cas3iterite  (SnO^)  which  has  the  same 
structure  as  rutile,  also  caused  a  significant  decrease  in  thermal  expansion 
anisotropy. 

Two  principal  methods  were  used  to  predict  the  thennal  expansion  of 
two-phase  ceramic  bodies  knowing  the  thermoelastic  properties  of  the  individual 
phases.  Earner's  method  was  applied  to  ceramic  bodies  for  the  first  time. 
Satisfactory  predictions  were  made  for  several  two-phase  systems. 


ABSTRACT 

Several  methods  were  used  to  predict  the  thermal  expansion  coefficients 
of  pure,  single-phase  ceramics.  The  predictions  were  extendod  to  cover  1700 
pure  phases.  Based  upon  these  calculations  a  number  of  phases,  predicted 
to  have  low  values  of  thermal  expansion  coefficient,  were  synthesized  and  the 
thermal  expansion  properties  of  these  materials  were  measured.  Several  of 
these  were  found  to  be  low-expansion  materials.  The  thermal  expansion 
properties  of  cubic  UPgO^  are  particularly  interesting.  This  material  expands 
with  increasing  temperature  up  to  about  U00°C.  Above  this  temperature  the 
crystal  contracts,  returning  to  its  room  temperature  length  at  about  1000°C. 
These  properties  are  believed  to  be  unique. 

Variation  of  the  thermal  expansion  anisotropy  by  addition  of  solid 
solution  atoms  to  several  base  crystals  was  attempted.  Addition  of  vanadium 
to  rutile  (TiO^)  resulted  in  a  significant  decrease  in  thermal  expansion 
anisotropy.  Vanadium  additions  to  cassiterite  (SnO^)  which  has  the  same 
structure  as  rutile,  also  caused  a  significant  decrease  in  thermal  expansion 
anisotropy. 

Two  principal  methods  were  used  to  predict  the  thennal  expansion  of 
two-phase  ceramic  bodies  knowing  the  therraoela3tic  properties  of  the  individual 
phases.  Kerner’s  method  was  applied  to  ceramic  bodies  for  the  first  time. 
Satisfactory  predictions  were  made  for  several  two-phase  systems. 


FOREWORD 


This  final  technical  report  covers  research  performed  on  a  program 
sponsored  by  the  Department  of  the  Navy,  Bureau  of  Weapons,  under  Contract 
NOrd-l8hl9«  The  report  includes  significant  sections  from  reports  previously 
submitted  for  the  period  June,  19J>8  to  September,  i960  and  complete  coverage 
of  research  performed  during  the  period  September,  I960  to  July,  1962. 

This  research  was  performed  under  the  general  technical  direction  of 
Mr.  S.  J.  Matesky  of  the  Bureau  of  Weapons. 
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I.  INTRODUCTION 

The  thermal  expansion  of  solids  is  one  of  several  phenomena  included 
in  the  general  concept  of  the  equation  of  state  of  solids.  Several  aspects 
of  the  thermal  expansion  of  ceramic  materials  were  investigated  in  this 
program.  Among  these  were  the  followings  •••  • 

1.  Prediction  of  the  thermal  expansion  properties  of  pure 
ceranic  crystals.  Crude  prediction  methods  useful  in 
the  search  for  low  expansion  materials  were  investigated 
along  with  more  sophisticated  methods  useful  for  calculating 
expansion  coefficients. 

2.  Investigation  of  the  thermal  expansion  anisotropy  of  solid 
solutions. 

3.  Prediction  of  the  thermal  expansion  coefficients  of  two- 
phase  ceramic  bodies,  knowing  the  thermoelastic  properties 
of  the  individual  phases. 

In  recent  years  rapid  progress  has  been  made  toward  understanding 
thermal  expansion  phenomena  in  terms  of  atonic  vibration  mechanisms.  The 
objective  of  this  research  has  bear  to  indicate  ways  in  which  this  knowledge 
can  be  used  to  solve  practical  problems.  In  addition,  several  new  low- 
expansion  phases  were  found  and  previously  unavailable  thermal  expansion 
and  elastic  property  data  were  determined. 
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II.  PREDICTION  OP  THE  THERMAL  EXPANSION  COEFFICIENTS 
OF  PURE  CERAMIC  CRYSTALS 

JU  Introduction 

Several  methods  can  be  used  to  predict  the  thermal  expansion  coefficients 
of  pure  ceramic  crystals.  The  choice  of  a  method  depends  to  a  great  extent 
on  the  amount  of  available  information  and  the  accuracy  required  in  the 
prediction.  Since  both  of  these  factors  may  be  subject  to  wide  variations, 
it  is  desirable  to  have  several  different  prediction  methods.  The  methods 
investigated  in  this  program  are  listed  in  Table  I  together  with  listings  of 
the  information  required  for  prediction.  The  predictions  made  in  the  program 
and  the  evaluations  made  of  these  predictions  using  our  measurements  and  those 
of  other  investigators  are  described  in  the  following  paragraphs. 


The  relationship  between  the  openness  of  the  crystal  structure  and  the 
thermal  expansion  coefficient  has  been  noted  by  several  authors. Several 
measures  of  the  openness  of  crystals  can  be  calculated  with  little  prior 
knowledge  of  crystal  properties.  One  such  measure  is  the  openness  ratio  (  R  ). 


in  which  V FYi  is  the  volume  of  a  formula  weight  of  the  solid 
and  Vj  is  the  volume  of  a  formula  weight  of  ions  or  atoms 

The  openness  ratios  of  a  number  of  materials  having  known  values  of  thermal 
expansion  coefficient  were  calculated  as  part  of  the  CAL  sponsored  research 
program  which  preceded  this  contract."*  The  effects  of  various  assumptions 
concerning  the  degree  of  ionic  binding,  on  the  volume  of  the  atoms,  were 
determined.  Calculations  using  Pauling* s^  electronegativity  concept  to 
determine  the  ionic  radius,  with  the  additional  assumption  that  there  is 
a  linear  variation  of  ionic  radius  with  percent  ionic  binding,  seemed  to 
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TABLE  I 

METHODS  USED  FOR  PREDICTION  OF  THE  THERMAL 
EXPANSION  COEFFICIENTS  OF  CRYSTALS 


METHOD 

INFORMATION  REQUIRED 

1.  OPENNESS 

CHEMICAL  FORMULA, SPECIFIC  GRAVITY 

2.  OPENNESS  OF  HOMOLOQOUS  SERIES 

CHEMICAL  FORMULA, SPECIFIC  GRAVITY 

STRUCTURE  TYPE 

3.  ATOMIC  SPACING  AT  MINIMUM  FREE  ENERGY 
ASSUMPTION  OF  INDEPENDENT  VIBRATIONS 

CHEMICAL  FORMULA .STRUCTURE  TYPE 

INTERATOMIC  DISTANCE  (LIMITED  TO  IONIC 
CRYSTALS  VI TH  ONE  BOND  TYPE) 

H.  ATOMIC  SPACING  AT  MINIMUM  FREE  ENERGY 
USING  EVALD'S  METHOD:  VIBRATIONS  NOT 
INDEPENDENT 

CHEMICAL  FORMULA, STRUCTURE  TYPE 

INTERATOMIC  DISTANCE  (LIMITED  TO  IONIC 
CRYSTALS  VI TH  ONE  BOND  TYPE) 

b  PI-1273-M-12 


give  lie  beat  correlation  between  the  openness  ratio  and  the  thermal 

coefSdmt.  The  results  for  some  ceranic  phases  are  plotted  in  Figures  1 

and  2  iron  data  given  in  Tables  II  and  III.  At  high  values  of  the 

ratio  there  is  a  good  correlation.  At  lower  values  of  openness  ratio 

factors  become  important.  For  example,  the  phases  falling  farthest, 

curve  co  the  high  R  side  are  those  which  have  relatively  weak  bends 

cations  aid  isolated  silicate  groups  (wollastonite,  pseudo-wollaatesltBfc 

forsterite,  etc.).  Those  falling  farthest  from  the  curve  on  the  lo®?  M 

are  liases  such  as  corundum  (Al^O^)  which  have  especially  strong  heads*. 

Therefore,  it  seems  likely  that  this  approach  can  be  used  most  smtsssaSiSSg- 

when  applied  to  the  search  for  low-expansion  materials. 

33*  more  rigorous  theoretical  approaches  to  prediction  of  the  thesnaS. 

expansion  of  ceramic  phases  provide  some  confirmation  for  these  stateness**. 
2 

Smyth  ami  others  who  have  investigated  the  thermal  expansion  of  Egterials 

which  have  low  or  negative  thermal  expansion  coefficients  have  afct dba&afi. 

the  loar  expansion  of  some  materials  to  the  greater  freedom  fra"  transverse 

vibrathcra  in  materials  having  open  structures.  The  observed  changes  la; 

vibration  frequency  for  internal  oscillations  ( vib ratior.3 within  grasps)  and. 

7 

lattice  oscillations  described  by  Ganesan  for  the  case  of  sodium  ddnafo 
and  sodhsn  br ornate,  indicate  the  greater  importance  of  the  lattice 
oscillations,  over  that  of  internal  oscillations  in  determining  Hie  thermal 
expansion  of  crystals  with  group  formation.  The  importance  of  the  3trgaBgjttt 
of  bonds  is  recognizable  in  expressions  derived  by  Gruneisen  ,  Kontoronta?, 


Cartz^,  Borelius^",  Kumar^'^  and  others. 


The  openness  ratios  of  many  phases  have  been  calculated.  The  results 
of  about  11;00  of  these  calculations  were  presented  previously  The 
results  of  additional  calculations  are  presented  in  Appendix  A  of  this 
report.  Based  upon  the  results  of  the  first  series  of  these  calculations, 
a  number  of  phases  having  open  structures  were  selected  for  synthesis  ant 
thermal  expansion  measurement.  Ihe  results  are  summarized  in  Table  IT. 

As  a  reszlt  of  these  measurements,  several  new  low-expansion  phases  are 
available. 
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Table  III 

THE  OPENNESS  RATIO  AND  COEFFICIENT  OF  THERMAL  EXPANSION  OF  SOME  SILICATES 


M 


CHEMICAL 

FORMULA 

CONSTITUTION  AMO  METRO#  OPfNNtSS 

RATIO 

MM  mill  COCFFICinl 

Of  tMMM.  UPMIlOa 
•c-l  MIC  1 

ui2ornio2 

MHtC«*ST. ,  OOIMOM.IIO,  KUUITf 
IMElFCMMTtO  (MU  OF  OOKaltO 
tnwu) 

0.17 

M  «  IO*7 

»-» o*e 

USTIO(IOOO)- 

»!«{ 

POlVCRfST.,  Mt  CaSlAM,  MONOCUNlC 
IRICRFfiROMCnt 

o.u 

>1  i  10*’ 

10*100*0 

001(0(1070) 

JIO.AI^j. 

ISIO| 

SINGLE  CRYSTAL.  KtTt,  NIIARONAL 
INTERFEROMETER 

0.11 

10  >  I0*7 

10-210*0 

trucooofo) 

C-fl.ti02  . 

POLVCRYST. ,  PSEURO-NOUASTONlTt  MONO- 
ClIRtC,  HlATOMf  ICR 

0.50 

M  •  IO*7 

100*100*0 

iiiot(nro) 

CaO. SiO| 

POLVCRYST.,  PSEURO-NQLL  ASTON  ITT  MONO* 

clinic.  $uirmm 

•.so 

M  «  I0*7 

io«->oo*c 

01007(11711 

CaO.SiOj 

POLVCRYST.,  NHL  AS  ION  IT  f 
liTCRFCRONCTtt 

o.st 

no  ■  io'7 

10-100*0 

MOCOOiMOl) 

ZCaO.SIOj 

POLY  CRY  ST. ,  KaO.StO|  RNOMOIC, 
OllATOMCTtt 

0,U 

IM  •  IO*7 

100-100*0 

01007(1171) 

2Cat.  AljOj. 
SiO, 

CO.AIjOj. 

»*o2 

POLYCmr.,  TETRAOONU,  OCNUllTC 
OllATQMCTCR 

0. 64 

01  ■  M*7 

ioo-wo*c 

01007(1171) 

POLYCRVST..  TtiaiNlC.  ANORTHlTC 
ItTCRFCROMCTtR 

0.02 

*1  >  IO*7 

10-100*0 

0CUCKIO70) 

CaO.MfO* 

isio2 

poiYCinrST..  MMoaiNic.  oiopsioc 
OILATONCTtR 

0.59 

oo  »  io*7 

100-100*0 

01007(1171) 

CaO.MfO. 

II02 

POLVCRYST..  OTNORNRMHC,  MONTI- 
CCUITC,  DlLATRMf T(t 

O.SR 

III  ■  IO*7 

ioo-ioo*c 

01007(1171) 

2  CaO.MfO. 

2Si«2 

POLYCRYST..  TCTfAOMRL,  UCRMAIITI 
OILATONCTtR 

0-  so 

00  I  IO*7 

100-100*0 

■1007(1171) 

JCaO.MfO. 

2SIOJ 

POLYCRYST.,  MONOaiNIC,  MfRWtNtTt 
0ILAT9MCTER 

0.59 

10.  >  IO*7 

IM-K»*C 

01007(1171) 

IFiO.liOj 

POLVCRYST.,  ORTMORNOMOIC.  FAYALITC 
DILATOMCTER 

0.01 

0>  .  IO*7 

»-100*C 

■1017(1171) 

LiJ).  AIjO*. 
2SI0| 

POUCRYST.,  NEK  AGONAL  fUCRYPTlTf 
OILATOMCTtl 

O.Oi 

-0  i  IO*7 

2o-xw*e 

OIUUT  1  MM 
MTC  1 

lljO.OIjOj.OIIOj  « 
LI  jO.IIjOj.IIIOj 

F0l7Ci71T.,  MMOai.lC,  FfTALtT£(7) 
OIUTOMfTd 

0.00 

0  >  IO*7 

io-ioo*e 

■MNa(iiM. 

1100)  aOTt  1 

MfO.SJOj 

POtYCRTST.,  MfliOaiNlC,  CLINOINSTATITE 

0.  so 

0»  i  IO*7 

I00-300*C 

RICtV(U79) 

2MfO.SIOt 

POLVCRYST..  ORTMORMQMflC,  FORSTER! TC, 
0ILA10MCTER 

0.55 

01  ■  I**7 

100-30*5'  C 

RISRYfl  179 

2Mff.2AI203 

POLfCRTST.,  90S  COROICRI TC  INTERFER¬ 
OMETER 

0.43 

II  <  IO*7 

25-100*0 

ecuci(io7o) 

Na2O.AI2Oj. 

2Si02 

POLVCRYST.,  TRICLINIC,  CARNEOIEITE, 
OILATOMETER 

0.03 

■  iio  ■  io"7 

20-J00*c 
.nuanacl(i  ion) 

Na20.  CaO.  Si  02 

POLVCRYST.,  C8IIC,  OILATOMETER 

0.50 

150  >  IO*7 

20-300*C 

MMMEL(IION) 

Si  Oj 

VITREOUS  SILICA,  TEL  Dll  CROSCOPf 

0.72 

7  >  IO*7 

2S-300*C 

NNimM0RE(l265) 

ZrOj.JiOj 

POLVCRYST.,  TETRAGONAL,  ZIRCON, 

0.59 

it  ■  trT 

25-300 *C 
KNITTEMORCf  1265) 

NOTE  I 

THE  LITHIUM  ALUMINUM  SILICATES  ARE  Nl OUT  KaOwN  TO  HAVE  LOW  VAiilCS  Of 
THEWAL  EXPANSION  COEFFICIENT.  THEY  A  HE  ALSO  NOTED  FOt  ID*  ?**JiES  Of 
MODULUS  Of  RUPTURE  *■»»  ARE  THOUGHT  TO  BE  DUE  TO  INTERNAL  CRAWS  f OSMff 
AS  A  RESULT  Of  ANISOTROPIC  THERMAL  EXPANSION.  THE  X-RAY  DETERMINATION  Of 
THE  THERMAL  EXPANSION  COEf flCIEN T/5-EUCRTPTITE  IS  REPORTED  IN  THE  /OLtCrf- 
ING  REFERENCE. 

F.H.  CILLERY  AO  C.A.  SU5*»  ‘THERMAL  CONTRACTION  Of  yS-EUCRfPTI  T£ 
(UjO.AIjOj^SiOj)  Sf  X-RAY  AMO  OILATOMETER  METHODS®.  J.  AM.  CERAM. 

SOC.  M2  (*)  176-77.  (1969). 
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Table  IV 

THERMAL  EXPANSION  OF  PHASES  HAVING  OPEN  STRUCTURES 


FORMULA 

DESCRIPTION 

OPENNESS 

RATIO 

THERMAL  EXPANSION 
COEFFICIENT  TEMP. 
Oc-lxlO7  RANOE°C 

REFERENCE 

AljOj.AtjOg 

SAMPLE  DECOMPOSED;  X-RAY  EVIDENCE  OF 
RON-CRYSTALLINE  MATERIAL. 

0.69 

H  1  K vlIH Ln  LI  W. 

(1959) 

VIOrilMCI  CT  11 

CdO.CoOj 

SAMPLE  DECOMPOSED 

(1959) 

3CoO. PjOj 

X-RAY  ANAL.  INDICATED  SAMPLE  WAS  MAINLY 
3CoO.P2Oj 

0.78 

60.8 

20-700 

KIRCHNER  ET  AL 
(1959) 

Lij0.2Si  Oj 

X-RAY  ANAL.  INDICATED  Lij0.2SI02  PLUS 
LijO.SIOj  AND  QUARTZ 

0.63 

t05 

20-700 

KIRCHNER  ET  AL 
(1959) 

MnO.SiOj 

10  X-RAY  DIFFRACTION  ANAL. 

0.6* 

83.3 

20-900 

KIRCHNER  ET  AL 
(1959) 

ThOj.SiOj 

X-RAY  ANAL.  INDICATED  ONLY  TH02.SI02 

0.66 

61.8 

25-1015 

MERZ  ET  AL 
(I960) 

3SnO.PaOs 

X-RAY  ANAL.  INDICATED  MAINLY  SSnO.PjOj 
WITH  SMALL  AMOUNTS  OF  SnSO,  AND  Na^PjOj 

0.7« 

96.7 

20-700 

K 1 RCIWER  ET  AL 
(1969) 

38a0.AI203 

X-RAY  ANAL.  INDICATED  PRESENCE  OF  BaO. 
Al203  AND  BaO.  NO  PATTERN  WAS  AVAILABLE 
FOR  3Ba0.AI203 

0.E7 

62.5 

20-1032 

MERZ  ET  AL 
(I960) 

Sr0.2AI203 

X-RAY  ANAL.  INOICATEO  MAINLY  Sr0.2AI203 
SOME  FREE  SrO  ANO  AljOj  ALSO  PRESENT. 

0.S2 

67.0 

25-1052 

MERZ  ET  AL 
(I960) 

Sr0.AI203.2S:02 

X-RAY  PATTERN  SIMILAR  TO  BaO.AIjOj. 

2 S.  Oj  HOMOLOG 

0.62 

52.3 

26-1000 

MERZ  ET  AL 
(I960) 

CaO, Cu0.4S» Oj 

X-RAY  PATTERN  INDICATED  ONLY  CaO.CuO. 
ISiOj 

0.81* 

29.7 

39-1015 

MERZ  ET  AL 
(I960) 

SrO, CuO.ISiOj 

X-RAY  PATTERN  INDICATED  ONLY  SrO.CuO. 
1SI02 

— 

19.0 

26-1013 

MERZ  ET  AL 
(I960) 

BaO.CuO.NSiOj 

X-RAY  PATTERN  INDICATED  ONLY  BaO.CuO. 
USiOj 

36.6 

38-1007 

MERZ  ET  AL 
(I960) 

uo2.P2o5 

X-RAY  ANAL.  INDICATES  CUBIC  UP20;.  WITH 

A  POSSIBLE  TRACE  OF  (U0)2  Pj07 

0.72 

18. 8 
-32.3 
-  6.66 

20-400 

400-1200 

29-1200 

MERZ  ET  AL 
(I960) 

•THIS  COMPOUND  WAS  SELECTED  FOR  MEASUREMENT  BASED  UPON  THE  OPEN  SHEET  SILICATE  STRUCTURE. 
THE  OPENNESS  RATIO  R  WAS  CALCULATED  AFTER  THE  MEASUREMENT  OF  THE  THERMAL  EXPANSION 
COEFFICIENT. 
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The  most  unusual  lew -expansion  phase  studied  In  this  program  Is  cubic 
UPjjO^.  Expansion  measurements  were  made  using  both  the  x-ray  and  di  la  tome  ter 
methods.  The  results  for  both  methods  are  presented  in  Figure  3.  The  material 


Hfm  3  EfERMAL  EXPANSION  OF  UPjOy 


expands  up  to  about  idJO°C  and  then  contracts,  returning  to  the  room 
temperature  dimensions  at  about  ICOO°C.  Agreement  between  the  two  methods 
of  measuring  the  thermal  expansion  is  good  up  to  about  UOO°C.  Above  bOO°C, 
the  dilatometer  gives  slightly  higher  results.  Between  900  and  1000°C  the 
curves  cross  again  so  that  at  !CG0°C  the  dilatometer  measurement  is  lower 
than  the  x-ray  measurement*  The  lade  of  good  agreement  at  1000°C  and  above 
may  be  due  to  the  presence  of  a  small  percentage  of  (UOj^P^O^.  The  presence 
of  this  material  would  affect  the  dilatometer  measurement  since  this  method 
measures  the  overall  expansion  of  both  phases,  whereas  the  x-ray  method  is 


sensitive  only  to  changes  in  the  cubic  UPgO^  P^a3a  anc*  n°t  other  biases. 
Comparison  of  the  x-ray  patterns  for  UPgO^  at  room  temperature  and  1000°C 
(the  temperature  at  which  the  unit  cell  of  UPgO^  ^as  a^jn03'*'  returned  to  its 
rocm  temperature  dimensions)  dicws  close  agreement  in  both  *d"  spacing  a  and 
intensities.  The  pattern  taken  at  1000°C  is  excellent  and  gives  no  indication 
of  structural  changes.  Only  two  weak  lines  not  characteristic  of  cubic 
UPgO^  were  observed.  This  type  of  expansion  followed  by  contraction  without 
a  phase  transformation  is  unique.  Curvature  of  the  thermal  expansion  curve 
in  the  other  direction  is  not  uncommon  (Si,  Ge,  InSb  and  vitreous  SiO^  at 
low  temperatures)  but  the  reverse  has  not  previously  been  observed. 

The  thermal  expansion  of  uranium  pyrophosphate  can  be  compared  with  the 
results  for  the  similar  titanium  and  zirconium  compounds  measured  by  ftarrison 
and  Hummel. TiP^O^  has  a  rather  typical  thermal  expansion  curve.  The 
material  expands  almost  linearly  from  rocm  temperature  to  1000°C  and  has  a 
thermal  expansion  coefficient  of  approximately  122  x  10”^/°C.  for  the 
temperature  range  from  50°C  to  1DOO°C.  These  crystals  are  weakly  birefringent. 
On  the  other  hand,  zirconium  pyrophosphate  (ZrP^O^)  expands  slightly  more 
than  the  titanium  compound  in  the  range  from  room  temperature  to  iiOO°C,  but 
above  this  temperature  the  thermal  expansion  coefficient  decreases  so  that 
little  expansion  is  observed.  These  crystals  are  isotropic.  Harrison  and 
Hummel^state  that  the  "normal*  zirconium  phosphate,  ZrP^O^,  has  a  reversible 
inversion  from  a  low  temperature  cubic  form  to  a  high  temperature  cubic  form 
at  300°C  and  dissociates  to  the  zirconyl  compound  (ZrO^PgC^  at  1$50°C. 

The  thermal  expansion  of  UPgO^  seems  to  be  related  to  that  of  ZrP^O^  in  that 
the  material  expands  to  temperatures  in  the  range  from  300  to  1|00°C  and  then 
there  is  a  change  in  the  shape  of  the  thermal  expansion  curve.  In  the  case 
of  the  thermal  expansion  coefficient  gradually  decreases  becoming 

negative  at  temperatures  in  the  range  from  U00-5>00°C  so  that  the  sample 

contracts  with  increasing  tenperature  in  the  range  from  l|00°C  to  1200°C. 

2 

Smyth  has  explained  the  negative  thermal  expansion  of  vitreous  silica 
and  the  high  temperature  forms  cf  quartz,  tridymite,  and  crystobalite  on 
the  basis  of  the  variation  with  tsiperature  of  the  transverse  vibrational 
frequency  of  the  oxygen  atoms  vibrating  between  two  silicon  atoms.  As  the 
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volune  decreases,  the  transverse  vibrational  frequency  of  the  oxygen  decreases 

dne  to  an  increase  in  the  high  order  repulsion  forces  thich  resist  restoration 

of  the  oxygen  to  its  original  position.  This  increase  in  the  high  order 

repulsion  forces  is  only  partially  offset  by  an  increase  in  the  electrostatic 

attraction  forces  which  tend  to  pull  the  oxygen  back  into  its  position. 

According  to  Smyth  there  is  always  the  possibility  of  negative  expansion  if 

one  or  more  of  the  lattice  frequencies  can  decrease  in  this  way  as  the  volume 

decreases.  Perhaps  the  negative  expansion  of  UPgO^  can  be  explained  on  this 

basis.  Vfyckoff1^  describes  the  structure  of  ZrPgO^  (which  is  the  sane  as 

the  structure  of  UP„0_)  as  made  up  of  an  undistorted  array  of  Zr*^  and 
I.-  £  l 

ions.  The  P^O^  ion  consists  of  a  pair  of  PO^  tetrahedra  which  share 
an  oxygen  ion  with  a  P-0  separation  of  l.$6  angstroms.  The  other  P-0 
separations  are  1.52  angstroms.  Each  Zr+^  ion  is  surrounded  by  six 
octahedral ly  distributed  oxygen  neighbors  at  a  distance  of  2.02  X.  Relatively 
large  holes  run  through  the  structure  bordering  the  oxygens  shared  by  two 
phosphorous  atoms.  Since  there  are  no  close  neighbors  in  the  transverse 
directions,  a  displaced  oxygen  atom  is  drawn  back  into  position  by  electro¬ 
static  attractions  of  its  pair  of  phosphorous  neighbors.  This  is  opposed 
somewhat  by  the  high  order  repulsion  of  the  same  phosphorous  atoms.  As  the 
lattice  shrinks  there  is  a  decrease  of  the  restoring  force  on  the  transversely 
displaced  oxygen  atoms  and  a  decrease  in  the  transverse  vibration  frequency. 
This  is  the  condition  required  for  negative  expansion. 

C.  Openness  of  Homologous  Series  Method 

Many  authors  including  Hummel^,  Austin"^,  Megaw^,  Rigby^and  Smyth ^ 
have  stressed  the  importance  of  considering  the  crystal  structure  type  in 
ary  attenpt  to  systematize  the  thermal  expansion  properties  of  materials. 

Since  structure  type  is  an  important  factor  in  determining  the  thermal 
expansion  properties,  it  is  likely  that  variations  among  the  phases  within 
one  structure  type  will  be  substantially  less  than  the  variations  observed 
in  groups  of  phases  having  several  different  structure  types.  Therefore,  the 
effect  of  variations  other  than  structure  type  can  be  studied  in  series  of 


compounds  having  the  same  structure  type.  For  this  reason,  an  attempt  was 
made  to  correlate  thermal  expansion  coefficients  and  openness  ratios  for 
homologous  series  of  compounds.  Chemical  homologs  are  compounds  having  the 
same  structure  and  similar  chemical  formulas j  for  example,  the  series  of 
cubic  compounds  TiPgO^,  ZrPgO^  and  UPgO^. 

The  thermal  expansion  properties  of  several  homologous  series  have  been 
21  22 

studied.  *  These  studies  were  handicapped  by  lack  of  data  for  many  phases. 

In  sane  cases  Msgaw’s  rule  (ce»/pi//g’t  in  which  jo  is  the  coordination  number, 

g  is  the  ionic  charge  and  A  is  a  constant)  was  approximately  obeyed 

(cubic  alkaline  earth  oxides,  oxides  with  the  CaF^  structure,  and  carbides 

with  the  NaCl  structure).  That  is,  within  homologous  series  the  thermal 

* 

expansion  coefficients  are  the  same.  In  other  cases,  the  presence  of 
displacive  transformations  or  other  processes  made  application  of  this  method 
of  uncertain  validity  (the  cubic  pyrophosphates  TiP^O^,  ZrPgO^  and  UPgO^). 

In  the  case  of  the  alkali  halides  there  is  a  reciprocal  relationship  between 
the  thermal  expansion  coefficients  and  the  openness  ratios  fbr  the  alkali 
metal  and  halogen  series.  This  is  shown  in’  Figure  U.  In  these  calculations 
100£  ionic  binding  was  assumed.  The  thermal  expansion  decreases  with 
increasing  openness  ratio  for  the  halogen  series  F”,  Cl“,  Br~  and  I~.  This 
relationship  is  also  observed  for  the  alkali  metal  series  Li+,  Na+,  K+,  etc. 
Exceptions  occur  for  the  fluoride  series  LiF,  NaF,  KF  and  CsF.  The  halides 
with  the  cesium  chloride  structure  (CsCl,  "CsBr,  C3l)  also  show  a  small 
decrease  in  thermal  expansion  coefficient  with  increasing  openness  ratio. 

The  phases  with  the  cesium  chloride  structure  have  high  values  of  thermal 
expansion  coefficient  relative  to  those  with  the  sodium  chloride  structure 
due  to  an  increase  in  coordination  number. 

The  cubic  rare  earth  sesquioxide3  represent  an  interesting  homologous 
series  because  variations  in  chemical  bonding  are  minimized  while  at  the 
same  time  there  are  substantial  variations  in  the  openness  ratios.  Available 
data  on  the  thermal  expansion  of  these  crystals  are  given  in  Table  V.  The 
openness  ratio  is  plotted  vs.  thermal  expansion  coefficient  for  two  temperature 
ranges  in  Figure  $.  The  thermal  expansion  curves  for  four  members  of  this 
series  measured  by  the  x-ray  method  are  given  in  Figures.  6-9.  It  is 
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Table  V 

THERMAL  EXPANSION  OF  CUBIC  RARE  EARTH  OXIOES 


0X1 OE 

SAMPLE 
DESCRIPTI  ON 

THERMAL  EXPANSION 
COEFFICf BIT  I07/*C 

TBIPERATURE 

REFERENCE 

Sc2°3 

CU8IC 

85 

0-999°C 

STECURA  A  CAMPBELL  ( I98l)12 

*2®3 

85 

0-993°C 

STECURA  A  CAMPBELL  ( I9«l)** 

88 

20-  IOOO°C 

WARSHAW  A  ROY  (196 1)60 

s.2o3 

• 

CUBIC 

88 

0-950°C 

STECURA  A  CAMPBELL  (1901  )M 

Eu203 

CUBIC 

78 

0-IOII°C 

STECURA  A  CAMPBELL  (I96I)*2 

ed2o3 

CUBIC 

80 

0-l062°C 

STECURA  A  CAMPBELL  (  I96I)62 

OyjO, 

CUBIC 

82.9 

25-l200°C 

BROUN  A  KIRCHNER  (1981)®' 

CUBIC 

82 

0-985°C 

STECURA  A  CAMPBELL  (I9«l)*2 

83 

30'-8408C 

PLOETC,  XRY3TYNIAX  A  DOMES 

( I957)7 1 

Ho203 

CUBIC 

. 

82 

0-969°C 

STECURA  A  CAMPBELL  (I98I)®2 

^r?®3 

CUBIC 

82.1 

25-  I200°C 

BROUN  A  KIRCHNER  (1981)®' 

CUBIC 

79 

0-I038°C 

STECURA  A  CAMPBELL  (I98I)®2 

t"2®3 

CUBIC 

82 

0-  I03I°C 

STECURA  A  CAMPBELL  (I98I)®2 

Yb,0, 

CUBIC 

82.1 

25-l200°C 

BROWN  A  KIRCHNER  (l96l)0i 

CUBIC 

81 

0-l02l°C 

STECURA  A  CAMPBELL  (I96I)92 

lu2®3 

CUBIC 

78 

0-l015°C 

STECURA  A  CAMPBELL  ( I9« I)82 
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Figure  5 


THERMAL  EXPANSION  COEFFICIENT  VS.  OPENNESS  RATIO 
OF  CUBIC  RARE  EARTH  OXIDES 
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UNEAR  EXPANSION  -  PERCENT 


TEMPERATURE  °C 

Figure  6  LINEAR  THERMAL  EXPANSION  OF  5c,  Ow 
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Figure  8  LINEAR  THERMAL  EXPANSION  OF  Brz  ^ 


evident  that  there  is  little  variation  in  thermal  expansion  coefficient 
among  these  phases.  /The  thermal  expansion  coefficient  of  SCgO^  increases 
more  rapidly  with  increasing  temperature,  than  that  of  the  other  phases. 

The  conclusion  to  be  drawn  from  these  observations  is  that  the  openness 
ratio  does  not  provide  a  method  for  predicting  the  thermal  expansion  within 
homologous  series  except  in  the  case  of  series  like  the  alkali  halides  in 
which  there  are  very  large  variations  in  ionic  radius,  ionic  polarizability 
and  openness  ratio.  In  series  like  the  cihic  rare  earth  oxides,  carbides 
with  the  sodium  chloride  structure,  alkaline  earth  oxides  with  the  sodium 
chloride  structure,  etc.,  little  variation  in  the  thennal  expansion 
coefficient  is  observed. 

D.  Prediction  of  the  Thermal  Expansion  Coefficient  Based  Upon  Atonic 

Spacing  at  Minimum  Free  Energy _ _ _ 

Our  study  of  various  methods  of  predicting  the  thermal  expansion  of 
materials  included  an  investigation  of  the  atomic  spacing  at  minimum  free 
energy  utilizing  statistical  mechanics.  In  this  research  some  relationships 
which  exist  between  the  atomic  structure  of  solids  and  thennal  expansion  are 
pointed  out.  This  approach  is  important  because  our  principal  approach  to 
the  prediction  of  the  thermal  expansion  using  the  openness  concept  does  not 
involve  direct  knowledge  of  the  atomic  structure.  Therefore,  it  seems 
likely  that  many  of  the  inconsistencies  in  the  predictions  will  have  to  be 
explained  in  terms  of  the  atomic  structure.  In  the  following  paragraphs 
are  presented  some  explanations  of  the  effect  of  crystal  structure  on  the 
thermal  expansion  and  a  summary  of  calculations  of  the  thennal  expansion 
values  for  some  alkali  halides  having  the  sodium  chloride  structure.  The 
complete  calculations  are  presented  by  Kirchner^  based  upon  Snyth's  research. 

The  purpose  of  this  work  is  to  show  how  the  structure  of  some  crystals 
and  glasses  can  be  related  in  a  qualitative  and  even  semi-quantitative  way 
to  the  expansion  characteristics  of  the  materials.  Two  different  methods 
were  used.  The  first  method  is  not  mathematically  rigorous  since  the 
vibrations  of  the  individual  atoms  are  considered  as  independent  of  each  other 
rather  than  treating  the  crystal  as  a  whole  and  computing  the  proper  modes  of 
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vibration.  The  latter  is  possible  but  extremely  time  consuming  while  the 

treatment  of  the  individual  atoms,  although  not  rigorously  correct,  can  be 

very  helpful  in  leading  to  a  better  understanding  of  the  relations  between 
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expansion  and  structure.  It  has  been  shown  1  that  a  mode  of  vibration, 
the  frequency  of  which  increases  with  the  expansion  of  the  lattice,  makes  a 
negative  contribution  to  expansion  and  vice-versa.  Some  of  the  ways  in 
which  this  principle  can  work  will  be  illustrated  below. 

1.  Independent  Lattice  Vibrations 

In  order  to  compute  the  frequency  of  vibration  of  an  atom  or  ion 
about  its  equilibrium  position,  it  is  necessary  to  know  the  forces  which  act 
on  it  when  it  has  been  displaced  by  a  known  amount  from  this  position. 

Since  the  atom  at  rest  is  in  a  position  in  which  all  the  forces  acting  on  it 
add  up  to  zero,  it  is  necessary  only  to  calculate  the  amount  by  which  the 
force  of  each  neighboring  atom  is  altered  when  the  atom  in  question  is 
displaced  by  a  known  amount.  In  order  to  allow  actual  calculation  of  some 
of  the  forces  involved,  the  interactions  are  treated  as  purely  ionic  with 
Coulomb  attractions  or  repulsions  modified  by  higher  power  repulsions.  The 
equations  for  calculating  the  restoring  forces  are  indicated  below.  These 
apply  when  the  magnitudes  of  the  displacements  are  small  as  compared  with 
the  interatomic  distances. 

If  an  ion  at  the  origin  has  a  neighbor  at  x, y,  z,  attracting  it 
with  a  force 

e-'r%  (II-l) 

in  which  ef  and  e#  are  the  charges  of  the  2  ions  and  r  is  the  distance 
between  than,  and  if  the  atom  at  the  origin  is  given  a  small  displacement 
to  u  ,  v  ,  w  ,  the  change  in  the  p(  -component  of  the  force  exerted  on  the 
displaced  ion  by  virtue  of  its  displacement  is 

(II-2) 
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in  which 


/'  =  X 


+■  Jf 


The  y  and  z  components  have  the  same  form.  If  the  ions  have  charges  of 
the  same  sign,  the  results  are  the  same  with  the  opposite  sign. 

If  the  ions  repel  each  other  with  a  force 


( II— 3) 

in  which  A  is  a  constant 

then  the  change  in  the  x  -component  of  the  force  acting  on  the  displaced  atom 
will  be 


A  J  ,  .  t/X  +VXV+WXZ 

{  U~(r+C  - p - 


) 


(II-4) 


The  y  and  x  components  have  the  same  form. 

The  effects  of  these  forces  can  be  illustrated  by  a  simple  example. 

In  Figure  10 an  ion  P  is  located  at  the  origin  between  a  pair  of  ion3  at 
{-a  ,  O ,  O  )  and  at  (a  ,0,0),  If  it  is  assumed  that  the  electrostatic 
attractions  between  P  and  Q  are  of  the  form  (1)  and  the  higher  order  repulsions 
are  of  the  form  (3)  then,  if  P  is  displaced  to  (  u,  Of  0  ),  the  force  acting 
on  P  by  virtue  of  its  displacement  will  be,  on  applying  (2)  and  (U) 


4+r  2ru+  *%'/*■  O. 


(XI-5) 


Figure  10  LOCATIONS  OF  THE  ATOM*  IN  A  SIMPLE  EXAMPLE  USED  FOR 
CALCULATIONS  OF  INTERATOMIC  FORCES 


If  the  ions  were  originally  at  their  equilibrium  separations,  the 
first  term  will  predominate  giving  a  positive  restoring  force.  As  a 
decreases,-  the  first  term  increases  more  rapidly  than  the  second  giving  a 
greater  restoring  force  per  unit  displacement  or  a  higher  frequency  of 
vibration.  In  this  case  the  restoring  force  has  come  from  the  high  order 
repulsions  diminished  somewhat  by  the  Coulomb  attractions. 

If,  on  the  other  hand,  the  ion  at  the  origin  i3  moved  to  (  o,  o  ) 
which  means  a  transverse  displacement,  the  force,  now  in  the  y  direction,  is 

2v  -  ^5*-  (n-6) 

If  the  original  separation  is  the  equilibrium  separation,  the  second 
term  will  predominate  and  there  will  be  a  positive  restoring  force.  As  a 
decreases  the  first  term  increases  more  rapidly  than  the  second,  diminishing 
the  restoring  force  per  unit  displacement  and  decreasing  the  frequency.  Here 
the  restoring  force  comes  from  the  Coulomb  attractions  diminished  by  the  high 
order  repulsions.  This  accounts  for  the  different  behavior  from  the  previous 
case.  A  vibration  of  this  kind  contributes  to  a  negative  expansion  and  one 
of  the  fbrmer  kind  to  a  positive  expansion. 

The  conditions  for  the  existence  of  a  vibration  of  the  type  leading 
to  a  negative  expansion  are  that  in  the  structure  there  must  be  some  atom3 
connected  only  to  a  pair  of  neighbors  about  180  degrees  apart  and  so  arranged 
as  to  have  considerable  freedom  of  transverse  vibration.  These  conditions 
are  met  in  several  of  the  forms  of  silica  and  other  materials  which  crystallize 
with  comparable  structures.  Vitreous  silica  shows  a  negative  expansion  at 
very  low  temperatures  and  very  low  expansion  at  room  temperature.  The  high 
temperature  form  of  cristobalite  shows  a  negative  expansion.  The  more  compact 
forms,  such  as  quartz,  where  the  more  efficient  packing  inhibits  the  freedom 
of  transverse' vibration,  show  positive  expansions. 

If  the  vibrations  of  individual  atoms  or  ions  are  considered  as 
taking  place  independently  of  each  other,  it  is  possible  to  develop  a  fairly 
simple  quantitative  treatment  for  crystals  of  cubic  symmetry.  The  examples 
considered  will  be  phases  having  the  sodium  chloride  structure.  Phases 
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having  this  structure  include  most  of  the  alkali  halides  and  crystals  of 
ceramic  interest  such  as  magnesium  oxide  and  other  alkaline  earth  oxides. 

In  the  sodium  chloride  structure,  a  cation  at  the  origin  is 
surrounded  by  six  anions  at  ( a,0,0),  (-<2,0,0),  (  O,  a-,  O)  t  {0,-a-,0)}  (  Oq#.) 
and  (OtOta)  where  a  is  the  closest  cation-anion  distance  and  is  one-half  of 
the  edge  of  the  unit  cube.  Equation  2  may  be  used  to  find  the  restoring 
force  on  the  ion  at  the  origin  when  it  is  given  a  small  displacement  to 
(u,V,W)  where  a  ,  V  and  W  are  small  compared  with  a  •  The  total 
effect  is  zero  and  in  the  same  way  it  can  be  shown  that  the  complete  group 
of  ions  at  any  one  distance  of  separation  from  any  one  of  the  cations  or  anions 
makes  a  zero  contribution  to  the  restoring  force  on  that  ion  -then  it  is 
displaced  in  ary  direction.  This,  of  course,  applies  only  to  that  part  of 
the  restoring  force  arising  from  the  Coulomb  interactions  and  applies  only 
in  the  case  of  surroundings  of  high  symmetry.  It  applies,  for  instance,  to 
the  vibrations  of  a  silicon  in  the  cristobalite  structure  but  not  to  an 
oxygen  ion  in  the  same  structure.  Therefore,  for  independent  vibrations  in 
crystals  having  the  sodium  chloride  structure,  the  restoring  force  on  any 
displaced  ion  can  be  computed  entirely  in  terms  of  the  high  order  repulsion 
forces. 

The  repulsion  forces  fall  off  very  rapidly  with  increasing  separation 
and  not  much  error  will  be  introduced  by  considering  only  interactions 
between  nearest  neighbors  and  second  nearest  neighbors.  If  the  force  of 
repulsion  between  any  ion  at  the  origin  and  ary  of  the  six  nearest  neighbors 
at  distance  a.  from  this,  is  of  the  form 


(H-7) 


then  the  X  -component  of  the  restoring  force  on  the  atom  at  the  origin 
arising  from  a  displacement  to  (  u ,  X,  W  )  can  be  written  using  the  equation 
ii.  Adding  the  six  equations  and  simplifying  gives 


(n-8) 
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If  77  is  greater  than  2  this  force  is  always  negative  for  a 
positive  displacement  u  indicating  that  it  really  is  of  the  nature  of  a 
restoring  force.  The  x  -component  of  the  restoring  force  depends  only  on  u 
which  means  that  the  vibrations  in  the  X  ,  y  and  z  directions  are 
independent  of  each  other. 

If  the  force  of  repulsion  between  the  ion  at  the  origin  and  any  of 
the  twelve  ions  at  a  distance  a  is  of  the  form 

7*  (ii -9) 

the  X  -component  of  the  force  on  the  atom  at  the  origin  when  displaced  to 
(  u ,  y  ,  W  )  can  be  determined  using  the  equation  U,  Adding  the  twelve  equations 
and  simplifying  gives 


Wf 


(2-n)u 


(n-io) 


Some  reasonable  estimate  now  has  to  be  made  for  A  f  ,  A2  and  71 
for  the  different  alkali  halides.  The  method  of  treatment  implies  the 
following  model.  Each  ion  has  three  independent  modes  of  vibration,  all  with 
the  same  frequency.  Just  enough  interaction  between  the  different  modes  is 
assumed  to  allow  the  attainment  of  thermodynamic  equilibrium  so  that  the  usual 
laws  of  statistical  mechanics  can  be  employed.  This  allows  the  calculation 
and  use  of  the  free  energy  of  the  system  and  the  condition  that  this  free 
energy  of  the  system  is  a  minimum  determines  the  state  of  exparsion  or 
contraction  of  the  crystal  at  any  given  temperature. 

The  energy  per  mole  can  be  written  in  the  form 


(u-ii) 
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where  y>j  is  the  quant™  number  describing  the  state  of  vibration  of  the  ith 

mode  of  vibration  of  the  cation,  and  T?  is  the  frequency  of  vibration  of  the 

cation,  rtj  and  T?_  describe  the  jth  mode  of  the  anion,  and  £s  is  the 

potential  energy  of  the  lattice  at  rest.  ^  and  1?_  are  all  functions 

of  the  cation-anion  separation.  In  the  usual  method  of  attacking  such  a 
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problem,  the  partition  function  Q  is  set  up  and  this  has  the  form  , 


Q*  <? 


~KT  )  e 


2KT 


,JA/ 


: ± 

,£KT 


[i-e^  f  1  ,-e& 


(11-12) 


The  Helmholtz  free  energy  U  is  given  by 

U  =  -KT  /*  Q 

Hence 


U*  +■  3NKT  —p.  +b(t-e 


(II-13) 


( II— II4.) 


In- the  absence  of  stress  or  external  pressure,  the  cation-anion  separation 
will  be  such  as  to  make  U  a  minim™.  The  quantities 

©,'  4r  “d  ®-'¥ 

are  the  Einstein  characteristic  temperatures. 


l~s  has  the  form 
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and  ftay  be  expanded  about  T  Q  >  the  value  of  T  for  -which  £  is  a  minimum 


M  is  the  Madelung  constant  for  the  lattice  and  B  is  a  constant 


characteristic  of  the  lattice,  PQ  is  given  by 


Me 


n  Be 


r  *  r  X+r 

'  O  rO 


(11-16) 


If ’/'*•  ro  +  AT  ,  then  takes  the  following  form  using  the  usual  Taylor 

series  expansion  and  omitting  terms  containing  powers  of  AT  higher  than 
the  second. 


*s- 


-A 'Me* 


(n-17) 


Here  AT  is  assumed  to  be  the  amount  try  which  the  value  of  the  cation- 
anion  separation  at  any  temperature  varies  from  what  thi3  separation  would 
be  in  the  absence  of  thermal  motion.  The  quantity  A  T  is  determined  by- 
choosing  that  value  which  minimizes  U  so  that 


di-18) 


The  substituting  of  the  derivative  of  £s  with  respect  to  AT  and 
solving  for  AT  gives 


N Me2 (n  -l)dr 


(II-l 9 
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Taking  the  derivative  and  simplifying 


Ar  - 


-JR  rr 


NMe*(n-t) 


L  ,  L  d®-  1  f  d®. 

2.  dr  2  dr  *  dr  *  c®~  ;  dr  J 


(11-20) 


But 


or  =• 


/  dfar) 
rQ  dT 


(11-21) 


and  if  d ^ dr  and  d  ©- / Jr  can  be  considered  as  constants 

independent  of  the  state  of  expansion  and  therefore  of  the  temperature,  then 
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(11-22) 


where 


-/©_) 
{  T  J 


is  the  Einstein  Function  of 


® 


In  order  to  compute  the  qiantities  and  ©_  ,  it  is  necessary 

to  know  the  frequency  of  vibration  of  the  cation  and  the  anion.  This  can  be 
determined  by  differentiating  the  potential  energy  functions  suggested  by 
Pauling  ^  to  give  the  force 


r  g*  Sare  (r++r_) 

F  re  - 


Jl-1 


,  77 +•  t 


(H-23) 
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for  the  cation  and  anion  and  equations  of  similar  form  for  the  forces  between 
cation  and"  cation,  and  anion  and  anion.  Comparing  these  with  equations  II-7, 
8,  9  and  10,  it  can  be  seen  that  by  equating  II-7  and  the  second  term  of 
11-23,  we  can  solve  for  Af  .  In  a  similar  way  we  can  find  A2  .  Then,  by 
substituting  the  values  of  A,  and  A2  in  II-8  and  11-10  and  adding  the 
two  terns  we  have  an  expression  for  the  restoring  force  per  unit  displacement 
as  follows. 

For  the  cation  this  is 


artfa-j)  ,  Sr)(7)-t)8Qe*(*r+)y,~i  C+  ,TT 

- Twtz - t - -  rro  uwu; 


<7*>; 


and  for  the  anion 


-l)60es(r+-t-r_)7>~1 1  3 rt(r-i) B0c*(zr.)”  '  C _ 

r  *+£  r  (fFr)”**  ~  "  r”*e 


(11-25) 


The  frequencies  of  vibration  of  the  cation  and  anion  are  given  by 


*  .  JL  fZSO. 
+  In  y rn+r"** 
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2tt  jm.r 
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where  7v+  and  77j_  are  the  respective  masses  of  the  cation  and  anion. 
Since 


pis  w  (II'27) 


and 


_  -A  I  ~C+  77+2.  1 
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and  £_ 


-5, 


Equations  11-24  and  29  can  be  used  to  evaluate  C+ 
evaluated  using  equation  13-13  from  Pauling^.  In  this  equation  SQ 


can  be 
must 

be  chosen  so  that  F  (p)  has  a  value  of  unity  for  p  ■  0.79.  Here  p 
is  the  radius  ratio.  When  this  is  solved.  It  gives  for  30  the  value  0.029 
for  77  “9.  This  choice  of  &0  allows  the  use  of  Pauling's  ionic  radii 
for  T+  and  7".  and  gives  the  proper  value  for  r  (the  cation-anion 
distance  for  all  of  the  alkali  halides). 

Combining  equations  11-22  and  30,  the  expansion  coefficient  can  be 


Using 


i  in  the 

form 

R(7)  +  2)ra  f  f 

NM&Z  (77-/J  l 

77 

=  9 

R 

-  8.314  x  10“7 

A7 

-  1.748 

e 

-  4.77  x  10"10 

(II-3D 


e.s.u. 

The  first  factor  which  is  common  to  all  the  alkali  halides  has  the 
value  711.  Using  this,  the  known  values  of  the  cation-anion  separations 
(  r0  )  and  the  values  of  the  Einstein  functions,  Table  VI  was  prepared 
showing  the  computed  expansions.  The  experimental  values  for  the  alkali 
halides  are  given  for  comparison.  These  are  the  "best"  values  from  the 
literature  as  given  by  Thielke.^l 
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Li  CALCULATED 

210 

320 

360 

EXPERIMENTAL* 

340 

438 

600 

Na  CALCULATED 

270 

380 

410 

EXPERIMENTAL 

360 

400 

415 

X  CALCULATED 

350 

430 

460 

EXPERIMENTAL 

370 

330 

400 

Rb  CALCULATED 

390 

460 

480 

EXPERIMENTAL 

- 

360 

375 

Ca  CALCULATED 

120 

- 

- 

EXPERIMENTAL 

320 

460 

470 

33 


PI-1273-M-12 


The  agreement  between  the  theoretical  and  experimental  values  seems 
to  be  quite  good  considering  the  assumptions  that  have  been  made.  It  may  be 
possible  to  improve  the  calculated  thermal  expansion  values  by  using 
experimental  values  for  the  exponent  of  the  repulsion  term  rather  than  the 
ninth  power  function  which  has  been  assumed. 

2.  Vibrations  of  the  lattice  As  a  Whole 

The  previous  computations  of  thermal  expansion  of  alkali  halide 
crystals  were  based  on  calculation  of  the  frequency  of  oscillation  of  the 
ions  under  the  assumption  that  the  ions  are  frse  to  vibrate  independently 
without  disturbing  neighboring  ions.  The  present  study  is  undertaken  to 
investigate  the  feasibility,  in  calculations  of  this  kind,  of  treating  somewhat 
more  rigorously  the  vibrations  of  the  lattice  a3  a  whole.  The  assumptions 
are  still  made  that  these  vibrations  are  simple  harmonic  vibrations  and  the 
polarization  of  individual  ions  is  completely  neglected.  Because  of  the 
great  volume  of  work  involved,  the  calculations  were  limited  to  six  alkali 
halides  which  are  lithium  fluoride,  lithium  chloride,  sodium  fluoride, 
sodium  bromide,  sodium  iodide  and  potassium  fluoride. 

3.  Method  of  Attack 

The  forces  between  the  ions  were  taken  to  be  those  used  by  Pauling 
in  his  calculation  of  the  interionic  separations  in  the  alkali  halides.  Hie 
high  order  repulsions  were  all  assumed  to  vary  so  that  the  potential  energy 
of  a  pair  of  ions  varies  inversely  as  the  ninth  power  of  the  separation  and 
the  constant  of  proportionality  was  taken  to  be  proportional  to  the  eighth 
power  of  the  sum  of  the  radii  ard  to  have  the  general  form  postulated  by 
Pauling.  This  completely  specifies  all  the  parameters  necessary  for  a 
description  of  the  dynamics  of  the  lattice. 

li.  Ionic  Forces 

The  role  of  the  Coulomb  forces  was  handled  by  the  method  of  Ewald 
which  gives  the  value  of  the  vector  at  any  ion  location  for  a  wave  described 
by  any  wave  number  vector  in  the  reciprocal  lattice.  The  method  was  modified 
to  give  directly  the  components  of  electric  field  instead  of  the  Hertz  vector. 
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The  high  order  repulsions  were  considered  by  calculating  the  forces 
arising  from  the  displacements  of  nearest  neighbors  and  second  nearest 
neighbors  only. 

The  frequency  of  a  wave  travelling  with  a  wave  number  vector  described 
by  a  point  in  the  reciprocal  lattice  was  computed  by  assuming  a  frequency 
for  this  wave  and  then  writing  down  the  equations  of  motion  for  each  of  the 
two  ions  in  each  of  the  x,  y  and  z  directions.  This  gave  six  equations  which, 
on  elimination  of  the  amplitudes,  gave  a  determinant  which  was  a  sixth  degree 
equation  in  the  square  of  the  frequency.  Most  of  these  deterrainantal  equations 
could  be  rather  easily  broken  down  into,  at  the  most,  second  degree  equations 
yielding  easy  solutions.  The  equations  were  all  solved  by  electronic  computer. 

S>.  Reciprocal  Lattice 

To  minimize  the  number  of  calculations,  sixty  four  points  uniformly 
distributed  in  one  unit  cell  of  the  reciprocal  lattice  were  taken  to  describe 
the  wave  numbers  of  travelling  waves.  Because  of  the  high  symmetry  of  the 
lattice,  this  involved  only  eight  different  calculations  each  yielding  six 
-values  for  the  frequency. 

6.  Expansion  Coefficient 

If  there  are  7)t  modes  of  vibration  with  frequency  modes 

with  frequency  ,  etc.  then  for  one  mole  of  the  alkali  halide  crystal 
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idle  re  A/  i3  Avogadro's  number. 
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A3  was  shown  by  Smyth  ,  the  coefficient  of  expansion  CC  of  the 
whole  crystal  will  be  given  by 
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where  R  is  the  gas  constant  in  ergs  per  mole  per  degree,  £3  is  the 
potential  energy  of  a  mole  of  the  crystal  when  the  ions  are  at  rest,  r~ 
is  the  cation-anion  separation  in  centimeters  and  £7X,)  is  the  Einstein 
function  of  X  . 

For  each  of  the  modes  considered,  the  frequency  P  was  computed 
d\> 

as  described,  above  and  was  found  by  increasing  r  by  one  percent, 

keeping  all  of  the  force  constants  the  same,  and  redetermining  each  frequency 
P  .  The  value  of  was  considered  to  be  given  by  • 

Equation  11-17  showed  that  Es  could  be  expressed  as  follows 
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where  Ar  is  the  amount  by  which  the  cation-anion  3  eparation  at  ary 
temperature  varies  from  the  equilibrium  value  it  would  have  in  the  absence 
of  thermal  motion.  From  this 
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If  the  cation-anion  separation  is  known,  as  it  is  for  each  alkali 
halide,  everything  is  known  to  determine  the  value  of  the  coefficient  of 
expansion  for  each  alkali  halide.  Table  VII  gives  the  values  of  the 
calculated  and  experimental  coefficients  of  expansion  at  300°K  for  each  of 
the  compounds  studied. 

7.  Discussion 

The  accuracy  of  the  calculated  expansion  coefficients  depends  upon 
a  number  of  assumptions  including: 

1.  Negligible  ionic  polarizability. 

2.  Simple  harmonic  vibrations. 

3.  Degree  of  interdependence  of  lattice  vibrations. 

4.  Choice  of  the  exponent  in  the  term  representing  the 
repulsion  forces. 
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Table  VII 

THERMAL  EXPANSION  COEFFICIENTS  FOR  THE  ALKALI  HALIDES 
AT  27°C  CALCULATED  AND  LITERATURE  VALUES 


COMPOUND 

CALCULATED  THERMAL  EXPANSION 

COEFFICIENT  BASED  UPON  EVALD'S  METHOD 
°C“I  *  I07 

EXPERIMENTAL 

DATA* 

LIE 

280 

340 

LI  Cl 

710 

436 

NaF 

340 

360 

NaBr 

630 

415 

Nat 

674 

470 

KF 

338 

370 

‘THIELKE'S  (1953)  BEST  VALUES 


PI-1273-M-12 


37 


The  choice  of  the  Ewald  method  for  allowing  some  interdependence  of  lattice 
vibrations,  requires  extensive  calculations  which,  in  'tarn,  limits  the 
number  of  results.  The  predicted  expansion  coefficients  are  of  the  correct 
order  of  magnitude  and  vary  in  accordance  with  the  experimental  values. 
However,  further  calculations  of  this  sort  for  the  alkali  halides  do  not  seem 
to  be  worthwhile. 

Variations  in  ionic  polarization  seem  to  account  for  many  of  the 

differences  between  the  calculated  and  the  experimental  expansion  coefficients 

(Table  VI).  The  influence  of  polarization  on  the  thermal  expansion  of  the 

25 

alkali  halides  has  been  studied  by  Weyl.  Using  Weyl's  approach,  we  note 
first  that  the  theory  gives  reasonably  correct  results  for  NaCl  and  NaBr. 
Substitution  of  Li  for  sodium  leads  to  calculated  results  that  are  too  small 
whereas  substitution  of  K  and  Rb  leads  to  calculated  results  that  are  too 
large.  The  principle  forces  determining  the  arrangement  of  the  ions  in  the 
structure  are  the  attractions  between  anions  and  cations  and  the  repulsions 
between  the  anions.  Neglecting  the  polarization  of  the  cations,  as  a  first 
approximation  we  note  that  polarization  of  anions  results  in  higher  electron 
densities  between  anions  and  cations.  At  the  same  time,  the  electron  density 
between  neighboring  anions  is  decreased.  Therefore,  the  effective  charge  is 
greater  for  coulomb  forces  between  anions  and  cations  and  less  for  coulomb 
forces  between  anion3.  As  the  temperature  increases,  the  electrons  become 
more  symmetrically  arranged  around  the  anions.  This  results  in  a  decreased 
effective  charge  for  the  anion-cation  attraction  forces  and  an  increased 
effective  charge  for  the  anion-anion  repulsion  forces.  Therefore,  the  lattice 
constants  increase  more  with  increasing  temperature,  than  would  be  expected 
without  polarization  of  the  anions  and  the  experimental  expansion  coefficient 
is  larger  than  the  calculated  expansion  coefficient.  Since  the  Li+  polarizes 
the  anion3  to  a  greater  degree  than  the  other  cations,  this  polarization 
effect  increases  the  experimental  values  for  the  lithium  halides  to  a  greater 
extent  than  the  other  compounds.  Within  the  group  of  lithium  halides  the 
polarization  effect  increases  with  increasing  anion  size  and  polarizability 
so  that  Lil  ha3  the  highest  expansion  coefficient  of  the  alkali  halides.  On 
the  other  hand  decreasing  the  polarizing  power  of  the  cation,  for  example 
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III.  THE  THERMAL  EXPANSION  OF  SOLID  SOLUTIONS 

A.  Introduction 

There  is  a  large  difference  in  the  thermal  expansion  coefficients  along 

26  27 

different  crystallographic  directions  in  many  anisotropic  materials.  } 

Single  phase  ceramic  bodies  prepared  from  materials  having  a  high  thermal 

expansion  anisotropy  tend  to  form  internal  cracks  upon  cooling.  The 

hysteresis  observed  in  the  thermal  expansion  curve  of  Ti0o  is  probably 

caused  by  these  internal  microcrack3  .  Other  polycrystalline  ceramic 

bodies  such  as  those  composed  of  MgO^TiO^,  AlgO^'TiOg  and  Li^O'Al^O^'PSiO^ 

(  / 3  -eucryptite)  have  low  thermal  expansion  coefficients  but  are  very  weak 

27 

because  of  a  high  expansion  anisotropy.  Reduction  of  the  thermal  expansion 
anisotropy  of  the  crystals  used  in  these  bodies  should  result  in  stronger, 
more  useful  bodies. 

The  effect  of  solid  solutiom  atoms  on  the  thermal  expansion  of  ceramic 

19 

materials  has  been  studied  in  a  few  cases.  Austin  has  suggested  on  the 
basis  of  literature  data,  that  extended  solid  solution  lowers  the  expansion 

•  29 

coefficient.  For  instance  Rigby,  Lovell  and  Green  showed  that  the  addition 

of  FeO-SiO^  to  CaO'SiO^  reduced  the  thermal  expansion.  On  continued  addition 

of  FeO*SiO„,  however,  the  thermal  expansion  coefficient  began  to  increase 
^  30 

somewhat.  Kozu  and  Ueda  have  found  in  their  study  of  plagioclase  that  the 
relative  amount  of  albite  and  anorthite  strongly  influences  the  magnitude 
of  the  thermal  expansion  anisotropy.  These  authors^  have  also  studied  the 
effect  of  impurities  on  the  thermal  expansion  of  diopside  minerals  from 
various  localities.  The  thermal  expansion  anisotropy  of  one  of  these  diopside 
crystals  tiiich  contained  a  relatively  large  amount  of  iron  was  found  to  be 
considerably  less  than  that  of  pure  diopside  (CaMgSigO^).  Ricker  and  Hummel^2 
have  reported  that  the  solid  solution  of  SiOg  in  TiC^  decreases  the  thermal 
expansion  coefficient. 

The  papers  cited  above  indicate  that  the  thermal  expansion  of  crystals 
can  be  varied  by  addition  of  solid  solution  atoms.  In  the  case  of  one 
naturally  occuring  mineral  sample,  variation  of  thermal  expansion  anisotropy 
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with  composition  was  noted.  The  objective  of  this  investigation  was  to 
study  the  variation  of  thermal  expansion  anisotropy  with  composition  for 
several  ceramic  solid  solutions.  Since  this  effect  has  not  been  studied  in 
series  of  prepared  solid  solutions,  it  was  decided  to  begin  a  study  of  the 
influence  on  the  thermal  expansion  anisotropy  of  the  replacement  of  Ti+^ 
in  rutile  with  ions  of  differing  size  and  charge  such  as  Zr+^  and  Kn+^. 

Ti0o  was  chosen  for  study  since  there  are  approximately  1*0  cations  which  are 
within  of  the  size  of  the  Ti  4  ion  so  that  substantial  solid  solution  is 
possible  with  many  ions.  Also,  rutile  possesses  the  relatively  simple 
tetragonal  structure  and,  therefore,  the  thermal  expansion  anisotropy  can  be 
determined  by  high-temperature  x-ray  diffraction  methods.  Later  in  the 
investigation,  samples  in  the  SnOg-V^O^  and  Al^O^-Cr^O^  systems  were 
investigated  to  determine  whether  these  solid  solutions  also  obey  the  rules 
that  appear  to  hold  for  rutile  solid  solutions.- 

B.  Experimental 

1,  Specimen  Preparation 

Except  for  Nb+^,^  Sn+^,^  Zr+^,^  and  Be+^,^  there  is  little 
information  in  the  literature  on  the  extent  of  solid  solution  in  TiOg.  Most 
of  the  additives  were  selected  for  the  present  study  on  the  basis  of  ion 
size  since  it  is  generally  considered  necessary  for  the  substituting  ion  to 
have  an  ionic  radius  within  about  of  the  substituted  ion  in  order  for 
extensive  substitutional  solid  solution  to  take  place.  Ti+^  has  an  ionic 
radius  of  0.68  X  so  that  the  radii  of  the  substituting  ions  should  lie  in 
the  rarge  of  0»f>8  -  0,78  X. 

Because  of  the  lack  of  information  on  the  extent  of  solid  solution 

in  TiOg,  10  mole  per  cent  of  the  oxide  additive  (based  on  the  cation)  was 

usually  used.  If  this  amount  of  substitution  takes  place  on  firing  the 

specimen,  a  relatively  large  change  in  lattice  constants  is  expected.  The 

samples  containing  0.5  and  1.0  mole  percent  additive  were  based  on  Johnson's 
37 

research.  It  should  be  noted  that  the  actual  extent  of  solid  solution  was 
not  known  in  most  cases. 
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A  standard  firing  tasperature  of  15C0°C  was  selected  on  the 
basis  of  the  NbgO^-TiO,,,  ZrC^-TiO,,  and  Be0-Ti02  phase  diagrams,  all  of 
which  show  considerable  solid  solution  at  this  temperature*  In  the  case 
of  the  additives  Li^CO^,  Sb203  ^2®$*  a  ^er!?era^ure  1200°C  was  used 

because  of  the  high  vapor  pressure  of  these  compounds  at  higher  temperatures. 

These  additives  were  all  reagent  grade  materials.  The  TiC>2  was 
Titanium  Pigment  Corporation’s  Titanox  TG.  The  analysis  of  this  material 
is  <  0.01*  CaO,  <  0,1%  W03,  <  0.008^  Mb,  <  0.CO7*  Sb^,  <  0.0i*£  SiOg,  0.C06£ 
Fe203>  10  ppm  V,  0.5  ppm  Mn,  5  PP”  Cu,  1;  ppm  Cr  and  10-20  ppm  Pb.  In  Table  VIII 
the  additives  are  given  together  with  the  relative  weights  of  TiOg  and 
additive,  and  the  firing  temperature  and  time. 

The  rutile  solid  solutions  were  prepared  by  ball-milling  the 
weighed  additive  and  Ti02  for  four  hours.  The  powder  was  passed  through  a 
U  mesh  standard  screen,  mixed  with  1$%  birder  solution  and  dried  at  110°C  for 
an  hour.  After  passing  the  dried  mix  through  a  16  mesh  screen,  it  wa3  pressed 
at  7000  psi  into  1"  discs  about  l/lt”  thick.  These  discs  were  fired  to 
temperature  and  then  cooled  to  room  temperature  in  less  than  30  minutes  to 
avoid  possible  precipitation.  The  x-ray  powder  specimens  were  prepaied  by 
crushing  the  discs  in  an  alumina  mortar  to  pass  a  325  mesh  screen.  A  sample 
of  pure  Ti02  for  use  as  a  standard  was  given  the  above  treatment.  The  SnOg- 
V20^  samples  were  prepared  as  described  above  using  a  firing  temperature  of 
1000°G  for  16  hours.  Also,  solid  solution  specimens  in  the  Cr203~Al203 
system  were  prepared  as  described  above  using  a  firing  temperature  of  l500°C 
for  17  hours.  Additional  experiments  were  performed  using  the  phases  VC>2, 

Ge02,  MgF2,  Mg0.2Ti02  and  AIgGj.TiOg.  These  latter  experiments  were 
unsuccessful  for  a  variety  of  reasons  such  as  decomposition  of  the  desired 
phases,  lack  of  stoichiometry  or  inability  to  get  the  compound  in  the  desired 
crystal  form. 


The  equipment  used  for  the  x-ray  thermal  expansion  measurements 
consisted  of  a  G,  E.  XRD-5  s pec tro goniometer  and  a  Tem-Pres  high  temperature 
x-ray  diffraction  furnace.  The  sample  holder  was  a  platinum  plate  with  a 
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Table  VIII 

DATA  ON  THE  PREPARATION  OF  RUTILE  SPECIMENS 


AODITIVE 

USED 

MOLE  % 
ADDITIVE 
(CATION  BASIS) 

WEIGHT  OF 
ADDITIVE  (g) 

WEIGHT  OF 
TI02  (g) 

FIRING 

TEMPERATURE 

°C 

MW 

H 

ZnO 

10 

11.3 

100 

■591 

mm 

u2co3 

10 

5.1* 

100 

MEM 1 

ID 

MgO 

10 

5.6 

100 

1500 

■9 

Cr2°3 

10 

10.57 

100 

1500 

Sum 

10 

20.25 

100 

1200 

■9 

HnOj 

10 

12.08 

100 

1600 

n 

V2°6 

10 

12.6* 

100 

1200 

■9 

W)3 

10 

33.2 

100 

1500 

KB 

W03 

1.0 

5.7 

19*. 3 

1500 

■9 

Mb205 

0.2 

*2 

2*0 

1*00 

IB 

Zr02 

10 

17.15 

too 

1500 

PbO 

10 

23.66 

76.3* 

1500 

17 

Y2°3 

0.5 

0.71 

1500 

17 

CaO 

0.5 

0.35 

■  . 

1500= 

17 

BeO 

1.0 

0.31 

1500 

17 

M0O3 

1.0 

3.58 

196.* 

1500 

17 

Si02 

30 

*8 

152 

1500 

15 

SnOj 

10 

16.75 

80 

1500 

17 

N 

0 

e 

« 

20 

37.7 

60 

1500 

17 

Sn02 

30 

6*  .7 

80 

1500 

17 

milled  depression  13  x  16  x  0.75  mm  deep.  Because  of  the  high  thermal 
conductivity  of  platinum,  the  maximum  difference  in  temperature  across  the 
sample  holder  was  found  to  be  2 °C  over  the  range  25-1000°C.  Using  a  precision 
potentiometer,  it  was  possible  to  maintain  the  temperature  within  +  0.5°C 
during  a  diffraction  peak  measurement. 

The  sample  holder  mounting  was  provided  with  rotational,  tilting  and 
translational  adjustments  so  that  the  specimen  surface  could  be  aligned  on, 
and  parallel  to  the  spectrometer  axis.  ■ 

Copper  radiation  (  A  *  1.5b0$0  X  for  )  was  used  for  all 
measurements.  The  x-ray  tube  was  operated  at  $0  KV  and  1 6  ma  with  a  0.0007" 
nickel  filter.  A  3°  beam  slit  and  0.2°  detector  slit  were  used  with  a 
scanning  rate  of  0.2°  (29)  per  minute  and  a  chart  speed  of  30  inches  per  hour. 
The  peak  positions  were  taken  as  the  midpoint  of  the  diffraction  curves  at 
one-half  peak  height.  Precision  in  measurement  of  the  peak  position  was 
within  +  0.01°(29).  To  attain  precision  and  accuracy  in  the  determination 
of  lattice  constants,  use  was  made  of  lines  in  the  back  reflection  region. 

These  were:  for  Ti02  (303),  29  =*  13k.7°  (£21)29  -  138.8°  and  (213)  29  -  119.7° 
for  Sn02  (323)  29  -  137.7°  and  (512)  29 'lltf.3*  for  A1 "g  (0,1;,  10)29  *  D£.2 
and  (2a6)29  =*  136.0°j  and  for  Cr^  (330)  29  -  137.5°,  and  (3,  2,  10)  -  Dtf. 8°. 

An  analysis  was  made  of  the  precision  and  accuracy  of  the  x-ray 

measurements.  At  2l*°C,  considering  diffraction  lines  for  130°  29,  or  above, 

random  errors  due  to  angular  measurement  and  specimen  centering  were  estimated 

to  be  within  +  0.000U  X.  This  was  checked  by  making  a  series  of  nine 

measurements  on  a  sample  of  Fisher  precipitated  silver  (99.99^  Ag)  usii^; 

the  (511)  and  (ij22)  lines.  All  of  the  measurements  fell  within  the  estimated 

precision.  The  accuracy  of  the  measurements  wa3  estimated  by  comparing  the 

extrapolated  ASTM  card  value  for  the  a  of  silver  at  2l*.0°C  (u.859  X)  with 

o 

the  observed  values.  Agreement  between  the  values  was  also  within  +  0.00 Oh  X. 
The  accuracy  was  therefore  within  the  limits  of  the  precision.  The  effect 
of  x-ray  absorption  on  accuracy  was  estimated  by  comparing  the  diffraction 
peaks  of  pure  silver  and  a  mixture  of  silver  and  graphite  in  the  ratio  l/3* 

The  mixture  had  an  aQ  value  of  0.0001  X  greater  than  pure  silver.  The 
absorption  difference  of  these  two  compositions  is  considered  to  be  much 


greater  than  those  between  the  oxide  solid  solutions  investigated  so  that 
errors  die  to  x-ray  absorption  are  considered  to  be  negligible.  At  high 
temperatures  the  data  are  less  precise  because  of  a  temperature  variation 
across  the  specimen  surface  of  +_  2.0°C.  This  is  equivalent  to  an  error  of 
■'+  0.00016  X.  Added  to  the  random  errors  0.0001;  X)  at  room  temperature, 
the  total  precision  is  +.  0.000 6  X  for  high  angle  lines  up  to  1000°C. 

Concerning  accuracy  at  high  temper ature,  systematic  errors  in 
temperature  measurement  and  sintering  of  the  sample  or  sample  movement  due 
to  expansion  of  the  furnace  parts  must  be  estimated.  Considering  the  latter, 
a  displacement  of  the  specimen  of  O.OOl"  was  found  to  change  aQ  of  rutile 
by  0.00006  X.  Assuming  the  maximum  movement  of  the  sample  surface  to  be  mils 
an  error  of  +  0.0003  X  i3  introduced  which  is  added  to  the  estimated  room 
temperature  accuracy  of  +_  0.000U  X  to  give  a  total  high  temperature  accuracy 

of  +  0.0007  X. 

To  obtain  an  overall  estimate  of  the  accuracy  of  the  measurements 

including  temperature,  the  a  of  silver  was  measured  to  600°C  and  compared 
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with  the  values  of  Hume-Rothery  and  Reynolds  (1938);  The  observed  values 
agreed  with  those  of  Hume-Rothery  within  +  0.0003  X. 

The  room  temperature  lattice  constants  were  measured  for  each 
rutile  solid  solution  specimen.  These  are  presented  in  Table  IX.  If  there 
was  a  significant  change  in  the  lattice  constants,  it  was  assumed  that  a 
considerable  amount  of  solid  solution  had  taken  place  and  x-ray  thermal 
expansion  measurements  were  made.  In  the  case  of  rutile,  the  peak  positions 
for  the  (303)  and  (521)  lines  were  measured  as  a  function  of  temperature  up 
to  1000°C.  These  data  were  then  used  to  obtain  the  variation  of  the  lattice 
constants  with  temperature. 

In  order  to  obtain  information  on  the  extent  of  solid  solution  of 

vanadium  in  TiO^,  a  series  of  samples  was  prepared  with  various  percentages 

of  vanadium.  As  can  be  seen  from  Figure  11,  the  a  lattice  constant  varies 

o 

linearly  with  composition  up  to  about  17-20  mole  %  vanadium.  The  results  are 
plotted  in  terms  of  the  compositions  originally  prepared  and  do  not  take  into 
account  changes  due  to  volatilization  of  A  chemical  analysis  of  the 

lli. 8  mole  percent  vanadium  sample  indicated  the  presence  of  16.1  mole  percent 
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Table  IX 

ROOM  TEMPERATURE  LATTICE  CONSTANTS  OF  RUTILE  SOLID  SOLUTIONS 


ADDITIVE 
HOLE  $ 


LATTICE  COHSTAKTS  A° 


PURE  TI02 

25 

9.5937 

2.9593 

.6992 

1  H+' 

29 

9.5996 

2.9587 

.6890 

10  Ll+I 

28 

9.5988 

2.9588 

.6990 

10  Pb+2 

22 

9.5933 

2.9598 

.6998 

9.2  Hb+S 

28 

9.8128 

2.9632 

.6829 

0.5  Y+3 

29 

9.5980 

2.959 

.6990 

0.5  C«+2 

26 

9.5933 

2.9589 

.6982 

1  Be+2 

27 

9.5938 

2.9689 

.6892 

10  Hg+2 

29 

9.5990 

2.9591 

.6981 

10  Zr»+2 

27 

8.5937 

2.9589 

.6991 

10  Sb+3 

27 

8.6957 

2.9829 

.6987 

10  V+5 

27 

8.5859 

2.9583 

.6951 

10  Cr+3 

28 

9.6989 

2.9580 

.6932 

10  Mn+2 

27 

9.5973 

2.9579 

.6939 

1 

28.5 

9.5998 

2.9586 

.6939 

10  H+6 

26 

9.5996 

2.9598 

.6882 

10  Zr*4 

21 

9.6180 

2.9899 

.6878 

I  Mo*4 

28 

8.5937 

2.9586 

.6990 

30  Si*4 

26.5 

9.5933 

2.9583 

.6951 

10  So*4 

26 

8.6099 

2.9817 

.6981 

20  So*4 

28 

9.6258 

3,0059 

.6998 

30  3n+4 

25 

9.6806 

3.0282 

.6525 

Figure  II  LATTICE  CONSTANT  a0  AT  VARIOUS  TEMPERATURES  vs.  COMPOSITION 
FOR  I200°C  FIRING  TEMPERATURE 
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vanadium  after  heating  to  form  the  solid  solution,  While  an  increase  in  the 
ratio  of  vanadium  to  titanium  is  unlikely  (an  error  in  chemical  analysis  is 
more  probable),  the  analysis  does  provide  assurance  that  a  substantial  fraction 
of  the  vanadium  does  remain  in  the  material  after  heating. 

C.  Discussion  of  Results 

In  discussing  the  results,  several  things  should  be  kept  in  mind. 

First,  except  for  the  case  of  vanadium  and  tin  in  TiO^,  the  extent  of  solid 

solution  in  the  titania  specimens  is  not  known.  Second,  the  valence  of  the 

additive  is  not  definitely  known  in  some  cases.  Some  work  ha3  indicated  that 

there  is  a  tendency  for  solid  solution  atom  to  assume  the  oxidation  state  of 

38 

the  host  lattice.  This  "valence  inductivity  *  effect  has  been  investigated  for 

the  case  of  low  concentrations  of  Mh  in  Al-O,,  where  the  Mh  at  first  assumes  a 

+•»  *  d 

+3  oxidation  state  in  imitation  of  A1  but  then  changes  to  a  +1;  state  at 

higher  concentrations.  Recent  electron  spin  resonance  experiments  reported 
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by  Gerritsen  and  Lewis  indicate  that  vanadium  replaces  titanium  and  is 

present  in  the  four  valent  state,  at  least  at  small  concentrations. 

The  question  of  valence  is  not  only  important  because  of  its  bearing 

on  ion  size,  but  also  because  electrical  neutrality  must  be  maintained  in 

the  lattice.  For  instance,  ions  with  a  valence  greater  than  four  would  require 

the  formation  of  Ti+^  ions.  Since  Ti+^  has  an  ionic  radius  of  0.7 6  8,  the 

overall  effect  might  be  similar  to  the  addition  of  an  ion  larger  than  Ti+^ 

(0.68  8).  Similarly,  the  addition  of  ions  with  valences  less  than  four 

requires  the  formation  of  oxygen  vacancies  (0- ^  radius  is  1.U0  8).  Mauer  and 

Bolz^  have  measured  the  thermal  expansion  anisotropy  of  an  oxygen  deficient 

sample  of  TiO_  with  the  composition  of  Ti01  up  to  a  temperature  of  IiOO°C. 

c  i  +-1 

In  such  a  sample  each  oxygen  vacancy  requires  the  formation  of  two  Ti  J  ions 
in  order  to  maintain  electrical  neutrality.  There  are,  then,  two  factors 
operating;  cations  of  larger  size  and  lower  charge  than  Ti+^  are  created  and 
many  anion  vacancies  are  introduced  into  the  lattice.  In  Table  X,  data  based 
on  Mauer  and  Bolz*  work  are  given  which  shew  that  the  Ti01  _  sample  has  a 
much  larger  expansion  anisotropy  (  A  CL  )  and  a  somewhat  decreased  volume 
expansion  coefficient  (GCV  )  and  average  expansion  coefficient  (C£av  )  than 
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Table  X 

THERMAL  EXPANSION  DATA  FOR  RUTILE  SOLID  SOLUTIONS 
(ROOM  TEMPERATURE  It)  |000°C) 


HOLE  % 
ADDITIVE 


LINEAR  EXPANSION 
COEFFICIENT 

1°  c 


EXPANSION 
AN  ISOTROPHY 
1°  C 


VOLUHE 
EXPANSION 
COEFFICIENT 
1°  C 


9  ROOM  TEMPERATURE  TO  400°C 
*  ROOM  TEMPERATURE  TO  800°C 


AVERAGE 
LINEAR 
EXPANSION 
COEFFICIENT 
1°  C 


PURE  TI02 

83.4  x  I0“7  ioz.7  x  |0"7 

24.3  x  IO"7 

274.5  x  IO"7 

91.5  x  IO"7 

io  v2os 

89.2  x  IO"7  94.8  x  IO"7 

5.4  x  IO"7 

273. C  IO'7 

91.0  x  IO"7 

10  Zr02 

77.7  x  IO"7  110.7  x  IO"7 

33.0  x  IO"7 

266.0  x  IO"7 

88.7  x  IO"7 

10  Sb203 

81.1  x  IO"7  103.0  x  I0“7 

21.9  x  IO"7 

265.2  x  IO"7 

88.4  x  IO"7 

9.2  Nb205 

82.5  x  IO"7  lll.l  x  IO"7 

28.6  x  IO"7 

276.1  x  IO"7 

92.0  x  IO"7 

io  *o3 

79.4  x  IO"7  105.7  x  IO"7 

26.3  x  IO"7 

264.5  x  IO"7 

88.2  x  I0“7 

30  SI02 

82.5  x  IO"7  104.8  x  IO"7 

22.3  x  IO"7 

269.8  x  IO"7 

89.9  x  IO"7 

10  Mn02 

77.6  x  IO"7  II 1.0  x  IO"7 

33.4  x  IO"7 

268.2  x  IO"7 

88.7  x  IO"7 

io  li2co3 

82.2  x  IO"7  108.0  x  IO"7 

25.8  x  IO"7 

272.4  x  IO"7 

90.8  x  IO"7 

10  SnOj* 

75.5  x  IO"7  104.0  x  IO"7 

28.4  x  IO"7 

248.6  x  IO"7 

85.0  x  IO"7 

20  Sn02* 

70.9  x  IO"7  93.0  x  IO"7 

22.1  x  IO"7 

234.8  x  IO"7 

78.2  x  IO"7 

30  Sf.Oj* 

58.4  x  IO"7  92.0  x  IO"7 

23.6  x  IO"7 

228.8  x  IO"7 

76.3  x  IO"7 

IOCr203 

75.1  x  IO"7  109.7  x  IO"7 

34.6  x  IO"7 

264.0  x  IO"7 

84.7  x  IO"7 

DATA  OF  MAUER  AND  BOLZ  ( 1 967) 

PURE  TiO,# 

78.6  x  IO"7  99.2  x  IO"7 

20.6  x  IO"7 

268.4  x  IO"7 

85.6  x  IO"7 

Ti0l.97* 

68.0  x  IO"7  104.5  x  IO"7 

36.5  x  IO"7 

240.5  x  IO"7 

80.0  x  TO"7 

PURE  TIOj# 

84.7  x  IO"7  105.0  x  IO"7 

20.3  x  IO"7 

274.4  x  IO"7 

91.4  x  IO"7 
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stoichiometric  TiO^.  In  Table  X  the  thermal  expansion  data  from  room 
temperature  to  1000 °G  are  listed  for  the  specimens  studied.  Those  ions 
which  produced  relatively  small  changes  in  the  lattice  constants*  Li,  Si, 

Sb  and  W,  also  produced  little  charge  in  the  thermal  expansion  anisotropy 
(  A  OC  ) ,  Except  for  Sn,  ions  which  increased  the  lattice  constants  (Zr,  Nb 
and  Mn)  also  increased  A0£  .  And  finally,  vanadium  which  decreased  the 
lattice  constants,  also  decreased  A&  to  a  remarkable  degree.  It  appears 
that  most  of  these  data  can  be  systematized  on  the  basis  of  ionic  radii}  if 
there  are  no  valence  changes  or  vacancies  formed,  the  substitution  of  an 
ion  of  larger  radius  increases  AOS ,  and  conversely,  an  ion  of  smaller  atomic 
radius  decreases  A  OS  *  According  to  this  interpretation,  Cr^+,  with  the  same 
ionic  radius  as  V+^  produces  a  larger  increase  in  ACS  because  of  the 
formation  of  anion  vacancies.  The  result  here  is  similar  to  that  obtained 
by  Mauer  and  Bolz  for  oxygen  deficient  rutile.  In  the  case  of  Sn+\  the 
ionic  radius  is  only  slightly  larger  than  that  of  Ti+^  and  there  is  almost 
complete  solid  solution  formation  between  SnO^  and  TiO^.^  Figures  12 
to  16  show  the  thermal  expansion  curves  obtained  for  pure  TiO^  and  for 
zirconium  and  vanadium  additions  to  TiO^. 

To  provide  additional  evidence  for  the  above  explanation  for 
changes  in  ACS,  further  work  was  done  in  the  Al^O^-C^O^  system.  These 
oxides  are  completely  miscible  in  all  proportions.  Unfortunately,  however, 
the  addition  of  only  a  small  amount  of  one  into  the  other  causes  a  large 
broadening  and  loss  in  intensity  of  the  diffraction  peaks.  Besides  pure 
AlgG^-Cr^O^,  only  two  solid  solution  specimens  were  measurable.  The  results 
are  indicated  in  Table  XI  and  Figure  17.  Since  the  Curie  temperature  of 
CrgO^  is  just  above  room  temperature  and  this  causes  a  discontinuity  in  the 
thermal  expansion  properties'*^  the  thermal  expansion  anisotropy  has  been 
calculated  for  the  temperature  range  U0C-1000°C.  The  observed  changes  of 
thermal  expansion  anisotropy  are  too  small  to  be  considered  significant. 

The  fact  that  the  sign  of  the  thermal  expansion  anisotropy  of  Al^O^  is  positive, 
while  that  of  Cr^O^  is  negative,  leads  to  the  conclusion  that  significant 
changes  in  thermal  expansion  anisotropy  will  be  found  for  larger  additions 
and  a  composition  with  zero  thermal  expansion  anisotropy  (defined  for  a 
particular  temperature  range)  will  be  found  in  the  intermediate  composition 
range. 
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Figure  12  LINEAR  THERMAL  EXPANSION  OF  PURE  RUTILE  (Ti  Oj) 
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LINEAR  EXPANSION  -  PERCENT 


LINEAR  EXPANSION  -  PERCENT 


LINEAR  EXPANSION  -  PERCENT 


PERCENT 


Figure  17  LINEAR  THERMAL  EXPANSION  OF  PURE  Cra  Oj 
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In  some  cases  polycrystalline  ceramics  made  from  crystals  having 
substantial  thermal  expansion  anisotropy,  show  thermal  expansion  hysteresis, 
lhis  effect  was  shown  for  rutile  by  Mauer  and  Bolz^0  and  Charvat  and  Xirgery.^ 
In  spite  of  this  complication,  the  thermal  expansion  coefficient  measured  ty 
Mauer  and  Bolz  by  the  dilatometer  method  was  90.3  x  10~^/°C  (0-1000°C)  which 
compares  favorably  with  the  average  thermal  expansion  coefficient  calculated 
from  the  x-ray  data. 

The  addition  of  vanadium  to  the  rutile  lattice  results  in  a 

progressive  decreases  in  thermal  expansion  anisotropy  with  increasing  vanadium 

content  up  to  10  mole  percent  vanadium.  At  10  mole  percent  vanadium,  the 

thermal  expansion  anisotropy  is  substantially  zero.  Continued  increase  in 

vanadium  content  results  in  a  reversal  of  the  anisotropy  so  that  at  low 

temperatures  the  percent  expansion  is  greater  in  the  direction  parallel  to 

a  than  in  the  c  direction.  An  increase  in  slope  of  the  thermal  expansion 
o  o 

curve  parallel  to  the  c  axis  at  high  temperatures  was  observed  in  some  cases. 
Little  change  in  the  volumetric  expansion  with  increasing  vanadium  content 
was  observed. 

According  to  Cronemeyer4-5  the  structure  of  TiOg  can  be  considered 
as  made  up  of  Ti-O-O-Ti-O-O-Ti  chains  oriented  normal  to  the  c  axis.  Therefore, 
replacement  of  titanium  by  vanadium  with  resulting  decrease  in  aQ  and  slight 
decrease  in  cq  can  be  thought  of  as  causing  a  decrease  in  the  length  of  the 
chains  but  having  little  effect  on  the  distance  between  them.  To  determine 
whether  or  not  vanadium  would  have  the  same  effect  in  other  crystals  of  the 
rutile  structure,  the  effect  of  vanadium  on  the  thermal  expansion  anisotropy 
of  cassiterite  (SnO^)  was  determined.  In  this  case,  too,  the  thermal 
expansion  anisotropy  was  reduced.  The  results  are  given  in  Table  XII  and 
Figures  18  and  19.  The  thermal  expansion  anisotropy  of  crystals  is  important, 
in  part,  because  of  the  effect  of  this  anisotropy  on  local  stresses,  porosity 
and  strength  of  polycrystalline  bodies  made  from  these  crystals.  At  high 
temperatures  during  sintering,  the  stresses  between  crystals  are  relieved. 
During  subsequent  cooling,  the  crystals  contract  unequally  causing  local 
stresses.  These  stresses  may  cause  local  failures  and  the  formation  of  pores 
that  are  large  in  two  dimensions.  This  change  in  microstructure  can  result 
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Figure  18  PERCENT  LINEAR  EXPANSION  OF  Sn02 
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in  low  strength  and  low  thermal  conductivity.  The  reduction  of  the  thermal 
expansion  anisotropy  of  rutile  and  other  crystals  may  make  it  possible  to 
prepare  ceramics  with  improved  properties. 
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17.  THE  THERMAL  EXPANSION  OF  TWO-PHASE  BODIES 

A. '  Introduction 

The  thermal  expansion  of  composite  materials  (including  multiple-phase 
ceramics)  can  be  predicted  from  the  themoelastic  properties  of  the  individual 
phases  providing  the  resulting  body  is  microscopically  continuous  and  the 
phases  do  not  enter  into  solid  solution.  The  degree  of  validity  of  the  theory 
is  dependent,  primarily,  upon  the  completeness  with  which  the  internal 
stresses  arising  from  thermoelastic  dissimilarities  are  described.  The  case 
in  which  the  composite  body  i3  homogeneous  and  macroscopically  isotropic  is 
of  immediate  interest  and  will  be  treated  in  this  report.  The  results 

22 

reported  here  are  an  outgrowth  of  research  reported  previously  by  Merz  et  al. 
In  this  earlier  research  the  inadequacies  of  available  prediction  methods 
became  apparent.  This  led  to  the  use  of  more  satisfactory  methods  which 
are  described  in  this  report. 

In  the  literature,  several  simplifying  assumptions  have  been  used  to 
permit  computation  of  the  linear  thermal  expansion  of  two-phase  ceramic 
bodies  from  existing  data.  These  assumptions  have  been  reviewed  and  a  mere 
complete  theory  for  multiple-phase  ceramic  bodies  is  presented. 

An  important  factor  in  improved  correlation  between  thermal  expansion 
theory  and  experiment  is  the  determination  of  the  intrinsic  elastic  properties 
of  each  ceramic  phase.  Gross  reduction  in  the  eLastic  moduli  due  to  finite 
porosity  resulted  in  attempts  to  produce  high-density  ceramic  bodies. 

Apparatus  was  assanbled  for  induction-heated  hot-pressing  to  obtain  high- 
density  spinel  and  mixtures  of  lyrex  glass  with  spinel.  The  Fyr ex-spinel 
system  has  been  studied  in  detail  both  theoretically  and  experimentally. 

B.  Relationships  Between  Elastic  Moduli  and  Stress  Wave  Velocities  In 


The  application  of  a  dynamic  stress  at  the  boundary  of  a  solid  will 
generate  a  stress  wave  (or  waves).  For  small  disturbances,  these  stress 
waves  will  propagate  away  from  the  source  with  a  velocity  dependent  solely 
upon  the  elastic  properties,  density,  and  geometry  of  the  medium.  Measure- 


ment  of  elastic  wave  velocities  in  appropriately  shaped  specimens  will 
provide  data  for  computation  of  the  elastic  moduli. 

Transmission  of  ultrasonic  pulses,  consisting  of  a  burst  of  radio 
frequency  (R.F.)  stress  waves,  provides  a  convenient  method  for  simultaneous 
measurement  of  stress  wave  velocity  and  attenuation  over  a  range  of 
frequencies.  High  frequency  compressional  waves  and  high  frequency  shear 
waves  in  ceramic  specimens  have  been  measured  using  the  electronic  test 
system  indicated  in  Figure  20.  The  stress  wave  velocities  (  C  )  in 
semi-infinite  media  are  given  by: 


where 


stress  wave  path  length 
stress  wave  transit  time 


and  by  Kolsky. 


where 


compressional  wave  velocity  Ct  ~ ^or  a>'k 


shear  wave  velocity. 


w 


Young's  modulus  (modulus  of  elasticity) 
shear  modulus  (modulus  of  rigidity) 
mass  density 

wavelength  of  stress  wave  in  medium 
cross  sectional  dimension  of  specimen 
Poisson's  ratio 


These  equations  along  with  the  relationships  between  the  isotropic  elastic 
constants,^  e.  g. 


2(/+yj 


3(1  ~2.y>) 
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Figure  20  8L0CX  DIAGRAM  OF  CAL  ULTSASOHIC  TEST  SYSTEM  FOR 
STRESS  WAVE  PULSE  TRANSMISSION  MEASUREMENTS 
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were  used  to 


derive  the  following 


t-y> 


Note  that  vhile  four  elastic  constants  may  be  of  interest,  only  two 
independent  elastic  constants  exist  for  isotropic  media,  hence  only  two 
different  stress  wave  velocities  need  be  measured. 

The  ultrasonic  pulse  method  also  permits  convenient  use  of  interferometric 
techniques  for  accurate  determination  of  temperature  coefficients  of  stress 
wave  velocity  and  attenuation.  By  simultaneous  excitation  of  two  similar 
specimens  maintained  at  different  temperatures,  small  changes  in  propagation 
characteristics  are  readily  observed  as  R.  F.  phase  changes.  Since  these 
specimens  may  be  excited  at  megacycle  frequencies,  detection  of  small  phase 
changes  constitutes  a  detection  of  velocity  changes  of  a  few  parts  per 
million  or  less. 


* 


Identical  relations  were  also  given  by  Birch. 
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The  attenuation  of  megacycle  stress  waves  in  ceramics  will  arise  from 
several  factors:  ' 

1.  Elastic  anisotropy  of  the  individual  grains  introduces 
scattering  losses. 

2.  Elastic  hysteresis  within  grains  introduces  irreversible 
absorption  losses. 

3.  Relaxation  effects  between  grains  due  to  grain  boundary 
motion  introduces  relaxation-absorption  losses. 

lr.  Thermal  conduction  between  regions  of  compression  and 
rarefaction  introduces  thermoelastic  losses. 

5.  The  presence  of  impurities  within  grains  introduces  an  atomic 
relaxation-absorption  loss. 

6.  The  presence  of  pores,  voids  and  microcracks  introduces 
another  form  of  acoustic  scattering  loss. 

All  together,  these  various  losses  are  called  internal  friction  losses  and 
are  intrinsic  to  the  material  (i.e.  internal  friction  losses  are  independent 
of  specimen  geometry).  Measurement  of  internal  friction  over  a  range  of 
temperatures  and  frequencies  provides  a  method  for  separating  the  several 
losses.  Grain  boundary  relaxation  generally  produces  an  internal  friction 
peak  at  slightly  elevated  temperatures  while  elastic  hysteresis  losses 
generally  increase  as  a  power  of  the  absolute  temperature.  Scattering  from 
microcracks,  however,  will  tend  to  decrease  with  temperature  since  the  size 
of  the  microcrack  decreases.  In  some  ceramics,  the  microcracks  seem  to 

■fc. 

completely  close  and  reheal,  hence  the  acoustic  scattering  loss  from  such 
scattering  centers  would  disappear  at  appropriately  high  temperatures.  If 
this  loss  mechanism  contributes  significantly  to  the  total  internal  friction 
at  lower  temperatures,  then  ultrasonics  provides  a  method  for  nondestructive 
determination  of  the  onset  of  microcrack3.  Such  information  would  be 
immediately  useful  in  computing  the  stresses  required  to  induce  microcracks. 

*  Private  discussion  with  Professor  W.  R.  Buessem,  Pennsylvania  State  University. 


Computation  of  the  thermal  stresses  'necessarily  involves  knowledge  of 
the  single  crystal  elastic  constants  for  each  of  the  ceramic  phases  alor^j 
with  their  temperature  dependence.  Such  data  are  generally  obtained  by 
ultrasonic  pulse  techniques  and  have  been  reported  for  a  few  ceramic  single¬ 
crystals,  e.g.  sapphire,  ruby,  magnesium  oxide,  and  magnetite  (see  Appendix  C). 
The  number  of  independent  elastic  constants  increases  from  3  in  a  cubic 
crystal  to  13  in  a  monodinic  crystal.  In  each  crystal  class,  the  several 
elastic  constants  may  be  used  to  compute  the  single-crystal  bulk  modulus. 

In  multiphase  ceramics,  the  bulk  modulus  will  be  a  function  of  the  porosity 
and  the  component  single-crystal  bulk  moduli  as  discussed  in  detail  in 
Appendix  C.  Since  each  of  these  properties  may  be  measured  independently, 
a  theory  of  elasticity  for  porous  ceramic  bodies  will  provide  the  basis 
for  computation  of  composite  elastic  moduli,  thermal  expansion  and  resulting 
thermal  stresses.  A  systematic  computation  of  thermal  stresses  in  multiphase 
ceramics  should  provide  the  insight  toward  the  design  of  high-strength 
ceramic  bodies. 

C.  Applications  of  Ultrasonic  Measurements 
1.  Stress  Wave  Velocity  Measurements 

The  ultrasonic  equipment  was  arranged  for  use  in  measurement  of  the 
onset  of  microcracks  in  polyphase  ceramic  bodies.  Because  of  the  desire  to 
measure  stress  wave  losses  introduced  by  the  smallest  possible  crack  dimension, 
megacycle  frequencies  were  chosen.  However,  at  megacycle  frequencies,  the 
acoustic  scattering  losses  due  to  both  elastic  anisotropy  of  the  grains  and 
the  presence  of  voids  or  pores  are  sufficiently  great  that  only  short  path 
lengths  through  the  specimens  may  be  used  successfully.  In  order  to  obtain 
sufficient  accuracy  in  the  determination  of  stress  wave  transit  time  in  bar 
type  specimens,  several  variations  of  ultrasonic  interferometer  techniques 
were  employed.  The  most  useful  technique  utilized  a  long  pulse-length  such 
that  the  pulse  reflected  twice  at  the  fused  silica-ceramic  specimen  interfaces 
as  indicated  on  Figure  21.  The  twice-reflected  pulse  produced  phase 
interference  with  the  directly  transmitted  pulse  introducing  combined  response 
nulls  and  peaks  as  the  carrier  frequency  was  varied.  Observation  of  the 
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frequericies  at  which  maximum  phase  cancellation  occurred,  permitted  computation 
of  the  phase  velocity.  The  expression  derived  for  stress  wave  velocity  is 
similar  to  that  given  by  McSkimin.^ 


2Lf 

n-h4r  + 


360 


where  fa  frequency  at  which  a  phase  interference  null  occurs 

phase  shift  at  the  specimen-buffer  block  interface 
in  degrees 

7 1 »  number  of  wavelengths  in  a  2L  path  length  of  the  specimen 


In  the  pulsed  carrier  measurement  system,  the  pulse  length  i3  always  chosen 
to  include  a  minimum  of  several  complete  R.  F.  cycles  in  order  to  minimize 
the  frequency  spectrum.  Thus  for  a  given  pulse  length,  a  minimum  frequency 
of  operation  exists  and  the  value  of  n  cannot  be  determined  by  counting 
from  zero  frequency.  However,  measurement  of  two  consecutive  null  frequencies 
provides  a  method  of  evaluating  both  71  and  C  .  For  specimens  whose 
geometry  is  selected  to  avoid  velocity  dispersion, 

2Lf  =  C  -  2L(f-bAf) 

n  +  -t-  2 ty'  y,  ±  ,  cr  4.2  ft 

■560  .  '  360 

and  77=  -£•  -  O.S  -  §-¥- 
Af  360 

where  A^  =  frequency  change  between  successive  nulls 

For  all  measurements  made  using  this  technique,  eleven  consecutive  null 
frequencies  were  measured  and  the  mean  frequency  change  between  nulls 
determined.  In  general,  the  stress  wave  velocity  was  accurate  to  three 
significant  figures  by  the  long  pulse  overlap  technique  titich  was  adequate 
for  use  in  thermal  expansion  computations. 

Shear  wave  velocities  were  measured  with  two  independent  methods: 

1.  Total  pulse  transit  time,  and  2.  The  long-pulse  overlap  method.  A 
summary  of  the  data  obtained  with  the  second  method  is  listed  in  Table  XIII. 
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A  length  of  tyrex  cane  was  cut  and  lapped  on  both  ends  to  obtain 
the  100^  tyrex  data.  For  this  specimen  it  was  convenient  to  measure  the 
null  frequencies  between  two  pulses*  the  difference  being  represented  by 
(/S +3P+/S  +  2Bi^2T-*‘2&z>  )  -  ( /S +/P+-3S +2&+2r-f-2Bs)  ~2P-2S 

where  T=  path  length  through  one  transducer 

S  =  path  length  through  one  silica  block 
P=  path  length  through  tyrex  specimen 

path  length  through  one  transducer/silica  bond 
path  length  through  one  transducer/specimen  bond 

It  can  be  shown  that  the  snail  effects  of  the  bonds  and  transducer  paths 
completely  cancel  out  with  the  interferometer  technique  jroviding  that  the 
reflection  interfaces  are  properly  identified.  As  a  check  on  the  accuracy 
of  these  measure'  ent3*  data  for  tyrex  from  two  other  investigations  were 
obtained  from  the  literature.  The  differences  between  these  three  sets  of 
shear  wave  velocity  measurements  were  about  1%.  The  long -pulse  null 
frequencies  were  generally  determined  to  only  1%  accuracy  (frequencies  were 
read  directly  from  the  signal  generator  dial*  higher  accuracy  may  be  obtained 
through  the  use  of  a  frequency  counter  when  necessary). 

The  total  pulse  transit  time  provides  a  convenient  check  on  the 
stress  wave  velocity,  but  the  accuracy  decreases  with  decreasing  specimen 
transit  time.  Typical  values  of  transit  times  are  given  in  Table  XIV  for 
the  shear  mode  transducer  assemblies  illustrated  in  Figure  20. 

If  the  total  transit  time  was  determined  with  +  3%  accuracy*  the 
ceramic  specimen  transit  time  would  be  known  to  only  +  k0%  for  the  values 
given  in  Table  XIV.  In  order  to  improve  the  attainable  accuracy*  a  fixed 
delay  line  having  a  delay  time  near  the  value  to  be  measured  was  used 
(Figure  20).  By  incorporating  a  delayed  trigger  circuit,  the  two  delayed 
pulses  could  be  observed,  simultaneously,  with  an  expanded  sweep.  With  this 
arrangement,  the  delay  difference  could  be  determined  with  +  3 %  accuracy. 
Thus,  with  the  fixed  delay  equal  to  the  time  delay  of  the  transducer-buffer 
block  assembly  (without  the  specimen),  the  delay  difference  becomes  equal  to 
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the  specimen  transit  time  and  nay  be  measured  directly  ■with  _+  yfi  accuracy. 
In  general,  all  three  methods  were  used  to  determine  shear  wave  velocities 
in  each  ceramic  specimen  in  order  to  minimize  operator  errors'. 


Table  XIV 

SHEAR  WAVE  TRAMSIT  TIMES  FOR  TEST  SPECIMEN  COMPONENTS 


SYMBOL 


2S  +  2Bt  +  28,  +  2T 


MEDIUM 

PATH 

LENOTN 

L 

INCHES 

SHEAR  WAVE 
VELOCITY 

10*  CM/ SEC 

SHEAR  WAVE 
TRANSIT  TIME 

MICROSECONDS 

FUSED  SILICA  BLOCK 

2.000 

3.760 

13.51 

INDIUM  SOLDER  BOND 

0.0001 

0.71 

0.003 

SALOL  BOND 

0.0001 

0.1 

0.006 

AC-CUT  QUARTZ 
TRANSDUCER 

0.00 125 

3.31 

0.0096 

CERAMIC  SPECIMEN 
(100*  SPINEL) 

0. 2505 

1.68 

1.36 

TRANSDUCER  f  BUFFER 
BLOCK  ASSEMBLY 

- 

- 

27.06 

TEST  SPECIMEN  ASSEMBLY 

- 

28.12 

Compressional  wave  measurements  were  severely  restricted  because  the 
transducer  assemblies  available  consisted  of  one  megacycle  barium  titanate 
ceramic  transducers.  The  frequency  response  of  these  transducers  was  poor 
at  the  pulsed  oscillator  carrier  frequencies  (5  megacycles  and  above).  Large 
response  ripples  due  to  reverberations  within  the  ceramic  transducers  precluded 
use  of  interferometer  techniques.  The  extended  rise  time  obtained  wi  th  these 
transducers  also  precluded  precise  determination  of  transit  times.  Preliminary 
measurements  of  ccmpresaional  wave  velocities  in  the  lyrex-spirel  series  were 
not  accurate  enough  to  provide  comparison  with  Kerner's  theory  for  Young's 
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modulus  and  bulk  modulus  of 'the  composite  bodies.  Fused  silica  blocks  and 
2 $  megacycle  x-cut  quartz  crystal  transducers  prepared  with  optical  polished 
surfaces  should  be  used. 


Elastic  Moduli  of  Pol 


irstalline  Ceramics 


Multiphase  ceramics  usually  have  average  thermal  expansion  coefficients 
which  lie  between  the  two  extreme  end  members.  Thi3  occurs  because  the 
internal  stresses  present  in  continuous  solid-phase  bodies  prevent  the 
high-thermal  expansion  phase  from  dilating  as  much  as  in  the  isolated  free- 
particle  state.  Of  course,  these  same  stresses  also  cause  the  low  expansion 
phase  to  dilate  more  than  in  the  isolated  state.  A  similar  effect  occurs  in 
simple  compression  of  composite  bodies  as  illustrated  by  the  two-phase  model 
in  Figure  22,  The  important  point  is  that  the  internal  stresses  are  present 
within  and  between  the  individual  grains.  Hence,  it  is  the  elastic  properties 
of  the  individual  grains  which  should  enter  into  the  internal  stress  calculations 
Or  restated,  it  is  the  single  crystal  thermoelastic  constants  which  d et ermine 
the  composite  themoelastic  properties.  Aside  from  the  anisotropy  considerations 
the  thermoelastic  properties  of  available  single-phase  ceramics  bodies  differ 
from  the  single-crystal  values  because  of  the  general  inability  to  produce 
fully  dense,  pure  single-phase  bodies.  Finite  porosity  causes  each  of  the 
three  principal  elastic  ihoduli  {S,K,G  )  to  decrease.  Lang^J  for  example, 
found  a  70%  decrease  in  Young’s  modulus  for  a  pure  alumina  body  with  28% 
porosity.  Because  of  these  gross  changes,  it  becomes  important  to  be  able 
to  predict  the  porosity  dependence.  In  fact,  porosity  accounts  for  most  of 
the  variability  found  between  reported  values  of  elastic  moduli  for  single¬ 
phase  ceramic  bodies. 


However,  a  comparison  of  static  and  ctynsrti0  moduli  requires  consideration 
of  two  additional  factors.  First,  the  measurement  of  elastic  moduli  under 
static  conditions  gives  values  for  the  isothermal  moduli  while  sonic  or 


ultrasonic  measurements  give  the  adiabatic  moduli.  The  difference  between 
these  two  moduli  is  dependent  upon  the  ratio  of  specific  heats  (  -p—-  )  and 

|  O  v>\/ 

is  generally  assumed  to  be  negligible.  Mason45  has  computed  the  elastic 


compliances  of  a-quartz  under  the  two  conditions  and  found  the  largest 
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a)  LINEAR  ARRAY  OF  GRAINS  WITH  UNIAXIAL  LOADING 


y 


b)  PLANE  ARRAY  OF  GRAINS  WITH  BIAXIAL  LOADING 


Figaro  22  DEFORMATION  OF  A  TWO-PHASE  MODa  SUBJECTED  TO 
COMPRESSIVE  LOADING 
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h.9 

difference  to  be  less  than  2%.  In  rutile,  the  differences  are  less  than 
0.3%,  For  isotropic  solids,  a  relation  given  by  Zener‘S  may  be  used  to 
evaluate  the  difference  between  the  adiabatic  Young's  modulus  (  V/x  , 
unrelaxed  modulus)  and  the  isothermal  or  static  modulus  (  Y$  ). 


Vu  cl2TVm 

vs 


(IV-1) 


where  Cv“  specific  heat  at  constant  volume 

fi  3  densiiy 

T  3  absolute  temperature 

OL  3  linear  thermal  expansion  coefficient 


For  anisotropic  solids,  a  relation  derived  by  Voigt  and  reported  by  Vick 
U9 

and  Hollander  may  be  used 


(IV-2) 


where  A 3  adiabatic  elastic  compliance  -  isothermal  elastic 
compliance 

C&i ,  OAfr  3  thermal  expansion  coefficients  in  the  i  and  k  directions 
Cp  3  specific  heat  at  constant  pressure 

12  2  n 

Using  the  following  values  for  MgO,  'i/x  =  3.03  x  10  dynes /cm  ,  7"  3  288  K, 

C  =  0.20  cal/g/°C,  <X  3  13$  x  10  V°C,  and  yo  3  3.58  g /ext?,  the  dynamic 
modulus  (  y/u  )  was  computed  to  be  0.$%  greater  than  the  static  modulus 
(  V&  ).  A  cursory  review  of  the  range  of  these  physical  constants  fbr 
other  singL  e-phase  ceramic  bodies  suggests  that  the  ratio  of  adiabatic  to 
isothermal  elastic  moduli  will  not  differ  from  unity  by  more  than  a  few 
percent. 

At  elevated  temperatures  the  test  specimen  may  deform  slightly  by 
creep  producing  an  apparent  lowering  of  the  modulus.  At  sonic  frequencies 
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and  above,  the  period  between  cycles  is  too  short  to  permit  stress  relaxation 
by  creep.  Hence,  sonic  measurements  provide  an  unrelaxed  modulus  which  is 
always  greater  than  the  static  or  relaxed  modulus.  Dynamic  elastic  moduli 
have  been  used  for  all  computations. 

Table  XV  is  a  compilation  of  the  elastic  properties  of  fully  dense 
oxide  ceramic  bodies  obtained  from  space  averaging  of  single -crystal  elastic 
constants  and  from  extrapolation  of  porous  body  data  to  zero  porosity.^ 


E. 


Theoretical  Thermal  Expansion  of  Multiple-Phase  Ceramic  Bodies 

Turner"^  considered  the  thermal  expansion  of  composite  bodies  by  assuming 
that  only  isostatic  stresses  occurred  within  each  phase  due  to  the 
dissimilarities  in  thermoelastic  properties.  Hence  only  the  modulus  of  volume 
elasticity,  or  bulk  modulus  (  K  ),  and  the  volume  coefficient  of  thermal 
expansion  ( <^v  )  were  required  to  define  the  composite  thermal  expansion,  i.e.. 


(X.\j  —  ^ 


(IV-3) 


where  *  volume  fraction  of  the  ith  phase 

a  volume  coefficient  of  expansion  of  the  ith  phase 
*  bulk  modulus  of  the  ith  phase 


For  isotropic  bodies,  the  linear  coefficient  of  thermal  expansion  (C6  ) 

i3  approximately  one-third  of  the  volume  coefficient  (<Xv  )*  so  that  equation 

(IV-3)  remains  unchanged  when  CL  is  used.  Previous  applications  of  Turner's 

22  53 

method  1  have  resorted  to  an  additional  simplifying  assumption  vhich 
introduces  considerable  error  in  many  systans,  i.e.  assuming  V  =Vt=  V^-Vi  t 
then  may  be  replaced  by  because  y=3K(/~2 v)  where  V  ■  Tourg's 
modulus  (modulus  of  elasticity),  V  »  Poisson's  ratio  and 


ZViVi 


(n-b) 


*  For  unrestrained  isotropic  materials,  / -HXAt  -(/-t-ccAtZ—O^JccAt) 
for  GC&t«f  ,  hence  CCy  —  3Cc,- 
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ELASTIC  PROPERTIES  OF  FULLY  DENSE  POLYCRYSTALLINE  CERAMICS  AT  ROOM  TEMPERATURE 


It  is  to  be  noted  that  for  the  entire  range  of  solid  materials,  Poisson's 
ratio  is  only  required  to  vary  from  zero  to  0«50  (the  latter  corresponds 
to  an  ideal  liquid).  Hence  a  small  numerical  difference  in  Poisson's  ratio 
may  reflect  a  large  difference  in  solid  state  properties.  The  agreement 
between  values  predicted  by  equation  (IV-U)  and  dilatometer  measurements  has 
generally  been  due  to  an  opportune  choice  of  values  for  the  elastic  constants. 
When  the  gross  effects  of  porosity  upon  elastic  moduli  are  ignored,  the 
resulting  apparent  disparity  of  elastic  moduli  data  has  tended  to  permit  a 
selection  of  those  values  which  gave  the  best  correlation  with  Turner's  method. 

An  examination  of  the  internal  stress  model  used  by  Turner  reveals  that 
shear  stresses  between  different  grains  and  between  similar  grains  with 
dissimilar  orientations  have  been  neglected.  Assuming  that  shear  stresses 
are  negligible  is  equivalent  to  stating  that  at  least  one  phase  behaves  as 
a  liquid.  Clearly,  such  a  model  is  not  adequate  for  use  with  ceramic  bodies. 

After  establishing  the  necessary  conditions  for  a  more  realistic 
internal-stress  model,  a  review  of  the  literature  revealed  that  Keraer'kad 
derived  expressions  for  the  thermoelastic  properties  of  composite  materials 
based  upon  such  a  model.  Kemer's  model  incorporates  both  the  volume  dilation 
and  the  grain  boundary  shear  stresses  which  arise  in  continuous-solid-phase 

•H* 

composite  bodies.  The  linear  thermal-expansion  coefficient  of  composite 
bodies  given  by  Kernel*  is 


CL 


Q ’-L  VI 


4-Gq 


-bC> 


(I7-S) 


*  Neither  Turner's  method  nor  Kerner's  method  is  applicable  to  discontinuous- 
solid-phase  bodies.  Thus,  the  presence  or  occurrence  of  microcracks 
decreases  the  validity  of  both  approaches. 
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where 


-53 

stands  for  a  summation  excluding  the  index  1 


Gy  Vi  V/ 

(7-5vt)Gt+(a-10vf)  15(1 -v,) 

J  /— »  i  w  \/. 


r  _r  {f-o  utA*  kjj  joy-xTi) 

'  T7-  OVA  .  V/ 

^(7-5-irjG,  +(Q-10is,)Gi  iSO-Vi) 


(17-6) 


Z^GKj 


_  *Lk£  - 

v3X/  +4G, 
3AC/  *  4G, 


(17-7) 


The  shear  modulus  (Go  )  and  the  bulk  modulus  (  KQ  )  refer  to  the  moduli  of 
the  composite  body,  Kerner's  method  requires  computation  of  these  two  moduli 
before  the  thermal  expansion  of  the  composite  bocfy  may  be  determined.  Since 
most  applications  of  thermal  expansion  data  also  require  data  on  elastic 
moduli,  these  intermediate  steps  serve  dual  purposes  (Figure  23).  Since  it 
is  generally  more  convenient  to  prepare  bodies  on  a  weight -fraction  basis, 
the  volume  fractions  may  be  replaced  by  -Q-  where  P{  is  the  wei^it  fraction 
of  the  ith  phase  and is  the  density  of  the  ith  phase.  Note  that  phase- 
density  does  not  effect  either  the  internal  stresses  or  the  composite  thermal 
expansion  except  by  defining  the  fractional  composition.  Obviously,  when 
the  densities  of  each  phase  are  equal  (e.g.,  in  the  MgO-HgO'A^O^  system), 
the  weight  fractions  become  equal  to  the  volume  fractions.  If  the  bulk 
moduli  of  each  phase  were  equal,  then  equations  (17-3)  and  (17-5)  would 
reduce  to  one  equation  corresponding  to  a  simple  rule  of  mixtures,  i.e.. 


ZcCiVi 


(17-8) 


*  Index  1  was  used  for  the  annihilated  suspending  fluid. 
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p.  Comparisons  Between  Theoretical  and  BnpiricaL  Thermal  Expansion 
in  Two-Riase  Materials  _ _ _ _ _  ■. 

1.  Aluminum-Silica  System 

Figure  2l*  illustrates  the  differences  between  predictions  by 

equations  (IV-3)  and  (IV-5)  for  the  aluminum-fused  silica  system.  Measurements 

<3 

of  thermal  expansion  for  two  different  compositions  as  reported  by  Kingery 
are  included.  The  difference  between  theory  and  experiment  may  be  due  to 
stress  induced  microcracks.  Microcracks  arise  when  internal  shear  stresses 
exceed  the  local  shear  strength  and  hence  the  occurrence  of  microcrack3  would 
cause  the  empirical  data  to  lie  between  the  values  predicted  try  equations  (IV-3) 
and  (IV-5)  as  was  found  for  the  aluminum-silica  system. 

Kingery  used  the  oversimplified  version  of  Turner's  method  (equation  IV-U) 
to  compute  thermal  expansion  coefficients  for  this  systan.  ‘  Since  Xoung's  moduli 
are  equal  for  these  two  materials,  the  assumption  of  equal  Poisson's  ratio  is 
equivalent  to  assuming  equal  bulk  moduli  and  equation  (IV-Ii)  reduces  to  a  simple 
rule  of  mixtures.  The  simple  rule  of  mixtures'  is  valid  for  a  compacted 
mixture  of  powders  in  which  no  binder  has  been  added,  no  sintering  accomplished 
and  hence  no  internal  stresses  developed.  This  model  does  not  seem  representative 
of  the  aluminum-silica  system  unless  all  aluminum-silica  grain  boundaries  have 
fractured.  Table  XVI  summarizes  the  differences  between  the  several  approaches 
to  theoretical  prediction  of  thermal  expansion  in  composite  bodies. 

2.  Magnesium  Oxide-Spinel  System 

In  the  magnesium  oxide-spinel  (MgO'Al^O^)  system,  the  difference 
between  predictions  by  equations  (IV-3)  and  (IV-5)  is  small  (Figure  25)®  In 
fact,  the  range  of  reported  values  for  the  linear  thermal-expansion  coefficient 
for  the  end  members  alone  constitutes  the  greatest  uncertainty  to  the  predicted 
values  for  the  two-phase  ceramic  bocft*.  Table  XVII  lists  the  values  used  for 
this  analysis. 

The  thermal  expansion  coefficient  of  several  other  compositions  in 
the  MgO-spinel  system  were  measured  with  a  quartz  tube  dilatometer.  However, 
the  reproducibility  of  the  measurements  was  unsatisfactory  (+  10 %),  These 
data  were  omitted  in  Figure  25  since  a  precision  of  +  1%  would  be  required  to 
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Figure  2¥  THERMAL  EXPANSION  IN  THE  ALUMINUM  -  SILICA  SYSTEM 
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LINEAR  THERMAL  EXPANS 


WEIGHT  PERCENT 


Figure  25  THERMAL  EXPANSION  IN  THE  MAGNESIA  -  SPINEL  SYSTEM 


validate  the  theory.  Even  with  the  several  obvious  improvements  in  the 
diiatoraeter,  it  is  clear  that  it  will  be  difficult  to  obtain  adequate 
precision  to  establish  the  relative  validity  of  equations  (17-3)  and  (17-$)  ^ 

in  the  MgO-spinel  system. 

The  concept  vhich  has  been  developed  suggests  two  perturbations  to 
the  thermal  expansion  coefficient:  1)  The  thermal  expansion  of  a  single¬ 
phase  polycrystalline  body  (having  the  same  chemical  composition  as  the 

single  crystal)  will  depend  upon  porosity  if  the  single  crystal  is  anisotropic 

* 

in  thermal  expansion  and  elastic  properties,  and  the  pore  phase  i3  continuous. 
This  dependence  should  arise  from  the  changes  in  the  internal  stress  pattern 
surrounding  each  grain  as  the  continuous-pore  phase  increases,  The  investigation 
of  porous  alumina  by  Coble  and  Kingery^^  was  made,  with  discontinuous  pore- 
phase  bodies  in  which  the  expected  change  in  internal  stress  would  be  small. 

They  concluded  that  porosity  had  no  effect  upon  the  thermal  expansion  in 
alumina.  2)  The  thermal  expansion  of  multiple-phase  bodies  will  depend  upon 
porosity  if  the  elastic  moduli  and  the  expansion  coefficient  are  dissimilar 
among  the  phases. 

These  perturbations  would  not  be  present  in -powder  x-pay  measurements. 
Dilatometer  equipment  available  for  this  study  did  not  provide  sufficient 
accuracy  to  conclusively  validate  these  predictions.  For  example,  several 
100%  spinel  bars  measured  during  this  program  produced  thermal  expansion 
curves  which  appeared  acceptable  but  which  were  not  reproducible.  The 
porosity-dependent  thermal  expansion  of  this  elastically  anisotropic  phase 
could  account  for  at  least  part  of  the  data  scatter,  if  large  scale  stresses 
induced  in  processing  result  in  local  stresses  between  grains  due  to  elastic 
anisotropy  (see  Section  I7-G). 

3.  Pyrex-Spinel  System 

Attempts  to  produce  bars  in  the  magnesia-spinel  systan  with  porosities 
of  five  percent  or  less  by  cold-pressing  and  sintering  were  not  successfxl. 

y  G 

*  The  ratios  of  and  based  on  data  given  in  Table  XV  provide 

a  measure  of  the  elastic  anisotropy.  Isotropic  crystals  would  have  a 
unity  ratio  for  each  modulus. 
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In  order  to  expedite  development  of  hot-pressing  equipment  to  produce  high- 
density  bodies  at  moderate  temperatures,  a  glass-ceramic  system  was  chosen. 
Pyrex  and  spinel  were  selected  as  representative  of  two-phase  systems  having 
grossly  dissimilar  thermoelastic  properties. 

■B* ** 

Glass  cane  was  crushed,  ball-milled,  3ieve-separated,  and  mixed 
with  minus  325-mesh  commercial  grade  spinel.  These  mixtures  were  hot-pressed 
in  an  induction-heated  graphite  mold  in  the  shape  of  4-inch  diameter  disks 
about  3/8  inches  thick.  Forming  data  for  the  three  disks  selected  are  listed 
in  Table  XVIII. 


Table  XVI II 

FORMING  CONDITIONS  FOR  HOT-PRESSED  PYREX-SPINEL  DISKS 


SAMPLE  NO. 

COMPOSITION 

BY  UEI GHT 

FORMINS 
PRESSURE 
(MUSE  READING) 

FORMING 

TEMPERATURE 

TIME  AT 
MAXIMUM 
TEMPERATURE 

as-i 

26*  PYREX 

6000  P3I 

mam 

60  MIN 

es-5 

60*  PYREX 

WOO  P3I 

mm 

30  MIN 

S3- 6 

76*  PYREX 

WOO  P3I 

30  MIN 

A  test  bar  0.2  x  0.5  x  4  inches  was  cut  from  the  central  portion  of  each 
disk.  Densities  of  the  rectangular  prism  specimens  were  measured  by  a 
displaced-mercury  weight-change  method  (Table  XIII). 

*  Pyrex  Code  7740,  Corning  Glass  Works. 

**  Particle  size  less  than  44  microns. 
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Thermal  expansion  of  the  Pyrex-spinel  series  of  bars  was  measured 
with  a  quartz-tube  dilatometer  over  the  temperature  range  from  20°  to 
•300°C  (Figure  26).  These  specimens  had  been  subjected  to  no  previous  thermal 
cycling  except  for  the  initial  cooling  after  hot-pressing.  Ko  further 
annealing  was  attempted  in  order  to  avoid  the  possibility  of  increasing  the 
number  of  microcracks.  Dilatometer  readings  were  taken  >diile  the  temperature 
increased  at  a  nearly  steady  rate  of  3°C  per  minute.  The  thermal  expansion 
was  found  to  be  a  nearly  linear  function  of  temperature  over  the  temperature 
range  of  100°  to  300°  C.  The  data  below  100°C  contained  apparatus-originated 
errors  and  was  excluded  from  further  consideration.  Contraction  measurements 
during  the  cooling  cycle  produced  a  hysteresis  effect  indicative  of  non¬ 
equilibrium  stress  relaxation  in  the  glassy  phase.  It  is  recognized  that 
future  measurements  on  glass-containing  bodies  will  necessitate  discrete 
measurements  at  stabilized  specimen  temperatures  even  when  the  test  temperatures 
are  well  below  the  glass  strain  point.  The  elastic  moduli  for  the  fyr ex- spinel 
system  as  computed  from  equations  (17-6)  and  (17-7)  are  shown  in  Figure  27. 

Figure  28  is  a  comparison  of  the  Pyrex-spinel  thermal  expansion 
measurements  with  the  predictions  obtained  by  equations  (17-3)  and  (17-5). 

The  superiority  of  Kerner's  method  is  clearly  demonstrated  in  this  system. 

Two  end  points  for  spinaL  are  shown.  The  upper  point  (<X  m  88  x  10~^/°C) 
corresponds  to  the  thermal  expansion  of  powder  specimens  measured  ty  x-ray 
diffraction  techniques.  The  lower  point  ((X-  ■  81  x  10~^/°C)  represents 
dilatometer  measurements  on  polycrystalline  bara.  The  difference  may  be  due 
to  elastic  anisotropy  in  the  spinel  crystal  (see  Section  IV-G).  The 
appropriate  value  to  be  used  in  prediction  of  thermal  expansion  in  Pyrex- 
spinel  mixtures  is  determined  by  consideration  of  the  origin  of  the  internal 
stresses.  These  stresses  arise  within  individual  grains  which  are  surrounded, 
on  the  average,  with  grains  having  the  elastic  properties  characteristic 
of  the  mixture.  Hence,  the  thermal  expansion  of  spinel  powder  (x-ray  technique) 
not  of  the  bulk  ceramic,  should  be  used  in  the  computations.  Improved 
correlation  would  be  expected  if  Kerner’s  modal  was  extended  to  include  elastic 
anisotropy.  Elastic  anisotropy  tends  to  lower  the  thermal  expansion  in 
porous  bulk  specimens  if  large  scale  stresses  induced  by  the  fabrication 
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LINEAR  THERMAL  EXPANSION  COEFFICIENT, 
«X  I07/°C,  AVERAGED  OVER  100°  TO  300°C 


Figure  28  THERMAL  EXPANSION  IN  THE  PYREX-SPINEL  SYSTEM 


PI-1273 


process,  are  present.  This  change  would  be  dependent  upon  the  relative 
amount  of  spinel  present  in  the  lyrex-spinel  system.  It  is  expected  that 
the  change  in  thermal  expansion  predictions  due  to  the  addition  of 
anisotropy  terms  would  be  similar  to  the  differences  between  the  two  curves 
labeled  "Kerner’s  Method"  in  Figure  27.  As  a  reasonable  approximation  (and 
one  which  greatly  simplifies  the  computations),  substitution  of  the  thermal 
.  expansion  coefficient  for  bulk  specimens  (dilatometer  data)  could  be  used  in 
equation  (17-5)  to  obtain  the  theoretical  thermal  expansion  in  a  two-phase 
system,  one  of  which  is  elastically  anisotropic.  This  simplification  would 
be  much  less  valid  if  more  than  one  phase  were  elastically  anisotropic. 

It  is  interesting  to  note  that  equation  (17-5)  predicts  a  thermal 
expansion  in  the  region  of  1$%  pyrex-85^  spinel  tfiich  is  slightly  higher 
(<X  *  88.5  x  10“V°C)  than  for  pure  spinel  (  cn  =  88.0  x  10-^/°C)  while 
-  equation  (17-3)  predicts  a  continuously  decreasing  thermal-expansion 
coefficient  with  increasing  Fyrex  content.  This  difference  serves  to 
emphasize  the  importance  of  shear  stresses  to  the  internal  stress  model  and 
to  computations  intended  to  define  composition  for  purposes  of  matching 
thermal  exparsion  coefficients. 

G.  Thermoelastic  Anisotropy 

In  single-phase  polycrystalline  bodies,  internal  stresses  may  arise 
from  either  thermal  anisotropy,  elastic  anisotropy,  structural  anisotropy 
or  combinations  thereof.  For  a  cubic  crystallite  imbedded  in  a  rigid, 
isotropic  bocy  with  a  matching  thermal  expansion  coefficient,  and  if  no 
restraints  are  present  at  the  sintering  temperature  and  no  temperature 
gradients  are  present  during  forming,  then  the  internal  stresses  on  the 
crystallite  reduce  to  zero.  However,  if  restraints  are  present  during 
forming  due  to  temperature,  pressure  or  compaction  gradients,  then  stresses 
may  arise  in  a  polycrystalline  medium  even  with  cubic  symmetry  due  to 
anisotropy  of  the  elastic  restraints. 
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For  a  spinel  crystal  of  the  type  (MgO*3.5  A1JD-)  reported  by  Verma,^ 

•ft  •  r 

the  anisotropy  ratio  is  2,1  to  1,  indicating  that  elastic  restraints  in 
a  poly  cry  s  tall  ine  bocfy  may  be  anisotropic  to  a  high  degree.  In  a  medium  of 
this  type,  the  anisotropy  superimposes  a  spatial  modulation  upon  the  internal 
stresses.  Konhomogeneo us  polycrystalline  bodies  will  generally  possess  regions 
of  significant  internal  stress  levels  (  e.  g. ,  cold-pressed  and  sintered  bodies 
in  which  uniaxial  pressure  was  employed  are  frequently  found  to  have  large- 
scale  laminar-type  cracks).  When  these  internal  stresses  are  significant, 
then  the  effects  of  elastic  anisotropy  will  be  significant  even  though  the 
thermal  expansion  coefficient  is  isotropic.  Hence,  for  precision  dilatometer 
measurements  upon  polycrystalline  bars,  it  is  not  sufficient  to  form  a 
specimen  in  a  manner  which  results  in  inhomogenities. 

In  general,  the  spatial  stress  function  for  a  crystallite  within  a 
polycrystalline  body  will  contain  terms  involving  the  direction  cosines 
up  to  the  fourth  power. ^  The  space  averaged  stress  obtained  by  integration 
of  the  cosi zb1  and  cosine^  functions  does  not  reduce  to  zero  for  elastically 
anisotropic  crystals.  Hence  the  thermal  expansion  of  single-phase 
nonhomogeneous  polycrystalline  bodies  may  be  expected  to  differ  from  the 
average  thermal  expansion  for  the  single  crystal.  Comparison  of  x-ray  and 
dilatometer  thermal  expansion  data  for  MgOAlgO^  spinel  indicates  that 
differences  are  often  found  (Table  XIX). 

H.  Elevated  Temperature  Effects 

Each  of  the  internal  stress  theories  requires  the  composite  body  to 
consist  of  a  continuous  solid-phase.  The  presence  of  microcracks,  from  any 
cause,  seriously  impairs  the  validity  of  these  several  predictions.  Most 

-*  Elastic  anisotropy  in  a  cubic  crystal  is  proportional  to  the  ratio  of 
elastic  compliance  constants  (  Sij  )  given  by  the  expression 

2.  fSa  -  Siz) 

where  a  unity  ratio  indicates  an  isotropic  body. 


(I 7-9) 


Table  XIX 

THERMAL  EXPANSION  OF  MgO.AI203  SPINEL 


TEMPERATURE  RANGE 


AVERAGE  LINEAR  THERMAL  EXPANSION  COEFFICIENT  X  I07/°C 


POWDER  SAMPLE,  X-RAY 
DIFFRACTION  TECHNIQUE 


BEALS 
ET  AL 
( 1967)  26 


ZIMMERMAN 
( I960)  65 


POLYCRYSTALUNE  BAR,  DILATOMETER 


COARSE  GRAIN 

FINE  GRAIN 

OTHER  SOURCES 
STUTZHAN 

ET  AL 

87 

(I960) 

WITTEMORE 

ET  AL  88 

THIS  STUDY 

62 

78 

87 

76 

90 

73 

8* 

90 

80 

- 

88 

86 

91 

- 

■  - 

- 

88 

98 

«a 

- 

25°  -  300°C 
26°  -  600°C 
26®  -  900° C 
26°  -  I000°C 
26®  -  I2D0°C 
26®  -  IROO°C 
26®  -  I500°C 


88.2 

87.8 


88.3 


83 

88 


90.5 


microcracks  arise  when  the  internal  stresses  exceed  the  local  strength  of 
the  body  during  cooling  from  forming  temperatures.  Hence,  prediction  of  the 
elevated-temperature  thermal  expansion  of  composite  bodies  (vhere  the 
occurrence  of  microcracks  would  be  le33)  should  provide  a  better  opportunity 
for  experimental  validation.  In  order  to  apply  equation  (IV -5)  to  elevated 
temperature  predictions,  the  temperature  dependence  of  both  the  elastic  moduli 
and  the  thermal  expansion  for  each  phase  must  be  known.  Wachtman,  et  al 
have  shown  that  Young ' s  modulus  for  several  oxide  ceramics  (e.g.,  sapphire, 
alumina,  magnesia,  and  thoria)  follow  an  exponential  temperature  dependence 
of  the  form 

Y*y„  -er.f  . 
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where  T  =  absolute  temperature  in  °K  ■ 

B,T0  are  empirical  constants 
'i fa  -  Young’s  modulus  at  0°X 

Based  upon  Wachtman's  data.  Young's  modulus  for  alumina  ceramics  would 
decrease  by  10^  from  2 5°C  to  1000°C.  Similar  data  for  other  ceramics  show 
a  more  rapid  decrease  with  increasing  temperature.  Smiley,  et  al3  have 
compiled  data  on  the  elastic  moduli  versus  temperature  for  spinel,  mullite 
and  several  other  ceramic  phases. 

A  few  observations  related  to  the  temperature -dependent  effects  can  be 
made.  Young's  modulus  for  spinel  decreases  with  temperature  more  rapidly 
than  for  MgO.  Assuming  both  the  shear  modulus  and  the  bulk  modulus  decrease 
in  a  similar  manner,  the  elevated-temperature  thermal-expansion  coefficient 
for  the  spinel-MgO  system  will  be  everywhere  higher  than  at  room  temperature. 
Both  increasing  (Table  XVII)  and  decreasing  wou^t*  contribute  to 

this  change.  An  increase  of  10  percent  was  obtained  by  hig)i -temperature 
dilatometer  measurements  for  the  5 0 %  MgO,  50 %  MgO’AlgO^  specimen  (Table  XVII). 

The  elastic  moduli  of  a  few  materials  increase  with  temperature  below 
800°C.  This  is  characteristic  of  high  silica  glasses,  for  example.  In  the 
aluminum-silica  and  Pyr ex-spinel  systems,  the  elastic  moduli  changes  may  be 
expected  to  produce  significant  differences  between  the  room  temperature  and 
elevated-temperature  thermal  expansion, 

I,  Conclusions 

Each  of  the  several  methods  used  to  predict  thermal  expansion  of 
composite  bodies  requires  data  on  the  elastic  properties  of  the  end  members. 

It  has  been  shown,  that  the  elastic  moduli  which  should  be  used  far  ceramics 
corresponds  to  the  values  for  the  fully-dense  polycrystalline  body  rather 
than  the  lower  values  associated  with  the  more  readily  available  porous 
ceramics.  The  greater  validity  of  Xerner's  method  over  Turner's  method  for 
prediction  of  thermal  expansion  for  multiple-phase  ceramics  has  been  clearly 
shown  for  three  different  two -phase  systems,  Tne  completeness  of  Xerner's 
model  for  internal  stresses  is  expected  to  reveal  the  superiority  of  Xerner's 
method  for  most  composite  bodies  (which  remain  continuous  solids). 


In  general,  the  thermal  expansion  coefficients  of  two-phase  bodies  will 
lie  between  the  end  members  for  all  ranges  of  composition  with  one  exception. 
When  the  thermoelastic  properties  are  grossly  dissimilar,  the  resulting 
thermal  expansion  may  lie  slightly  outside  the  bracketed  range  for  compositions 
which  are  predominantly  one  phase  (e.g.,  1%%  Pyrex-Q5%  spinel). 

Preliminary  calculations  for  polycrystalline  specimens  composed  of 
elastically  anisotropic  crystals  indicate  that  the  presence  of  structural 
inhomogenities  will  introduce  internal  stresses  which  will  tend  to  lower 
the  thermal  expansion  coefficient.  Under  these  conditions,  the  bulk  thermal 
expansion  for  polycrystalline  bodies  will  differ  from  the  single-crystal  or 
the  x-ray-measured  thermal  expansion  for  powder  specimens.  These  local 
stresses  may  arise  even  though  the  single-crystal  thermal  expansion  coefficient 
is  isotropic,  provided  that  large  scale  stresses  induced  during  fabrication 
are  present. 

The  elastic  properties  (e.  g. ,  K ,  G,  V,  zr  )  of  composite  bodies, 
including  multiple-phase  ceramics,  may  be  predicted  (with  small  error  for 
bodies  which  remain  continuous  solids)  by  a  method  developed  by  Kerner.  This 
study  appears  to  be  the  first  application  of  Kerner’ s  method  to  ceramic 
bodies. 


V,  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  Conclusions 

1.  The  relationship  between  the  openness  ratio  and  the  thermal 
expansion  coefficient  for  a  number  of  oxides  and  silicates,  was  investigated. 
In  these  groups  of  materials,  phases  having  high  values  of  the  openness  ratio 
usually  have  low  thermal  expansion  coefficients. 

2.  Several  phases  having  high  values  of  the  openness  ratio  were 
synthesized  and  the  thermal  expansion  coefficients  were  measured.  The 
following  phases  have  especially  low  thermal  expansion  coefficients: 

CaOCuO^iOg 
SrOCuOliSiOg 
Ba0*Cu0*ltSi02 

UP20^  (cubic) 

3.  Cubic  UP  0_  has  unique  thermal  expansion  properties.  The  crystal 
expands  with  increasing  temperature  up  to  about  1;00  C.  Above  this  temperature 
the  crystal  contracts,  returning  to  its  room  temperature  dimensions  at  about 
1000°C.  Since  the  peak  positions  and  intensities  of  the  x-ray  diffraction 
patterns  are  very  nearly  the  same  at  room  temperature  and  at  1000°C,  these 
expansion  properties  are  not  the  result  of  a  major  structural  change. 

U.  Little  variation  of  the  thermal  expansion  coefficient  with  openness 
ratio  is  observed  within  homologousseries.  Ihese  observations  confirm 
Megaw's  rule  / _,  _  Ap2)  . 

'  ?*  / 

3.  The  thermal  expansion  coefficients  of  the  alkali  halides  were 
calculated  based  upon  the  change  in  lattice  spacing  with  temperature  at 
minimum  free  energy  using  the  techniques,  of  statistical  mechanics.  These 
calculations  require  knowledge  of  the  atomic  structure,  ionic  charge  and 
atomic  spacings  at  one  temperature.  The  repulsion  exponent  was  assumed  to 
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have  one  value  for  the  entire  group.  Comparison  of  the  results  of  these 
calculations  with  experimental  values  gave  encouraging  results.  Improvements 
can  probably  be  made  by  taking  into  account  the  polarization  of  the  ions  and 
variation  of  the  repulsion  exponent. 

6.  The  thermal  expansion  anisotropy  of  crystals  can  be  varied  by 
addition  of  appropriate  solid  solution  atoms.  Addition  of  vanadium  atoms  to 
replace  titanium  in  rutile  (TiOg)  results  in  a  marked  decrease  in  die 
thermal  expansion  anisotropy.  At  an  addition  of  about  10%  vanadium  the 
thermal  expansion  anisotropy  reverses  so  that  the  axis  becomes  the 
direction  of  highest  thermal  expansion  coefficient.  The  effect  of  vanadium 
on  the  thermal  expansion  anisotropy  is  not  restricted  to  rutile  since  a 
similar  reduction  in  thermal  expansion  anisotropy  was  observed  when  vanadium 
was  added  to  cassiterite  (SnO^)  vhich  also  has  the  rutile  structure. 

7.  The  thermal  expansion  anisotropy  of  Cr^O^  is  opposite  in  sign  from 
that  of  corundum  (AlgO^).  Since  corundum  and  Cr,^  form  a  continuous  solid 
solution  series,  an  intermediate  composition  must  exist  for  vhich  the  thermal 
expansion  anisotropy  is  zero. 

8.  Kerner's  method  was  used  to  predict  the  thermal  expansion 
coefficients  of  two-phase  ceramic  bodies.  The  predictions  are  satisfactory 
and  the  superiority  of  Kerner's  method,  which  accounts  for  shear  stresses 
was  demonstrated. 

9.  In  cases  in  ihich  the  thermoelastic  properties  of  the  pure 
phases  are  grossly  dissimilar,  the  thermal  expansion  coefficients  of  some 
intermediate  compositions  may  be  slightly  outside  the  range  between  the 
end  members. 


APPENDIX  A 


"Openness"  Ratios  of  Some  Ceramic  Phases 

In  the  following  list  the  compounds  are  mainly  covalent  in  character 
of  chemical  binding.  Therefore,  the  calculation  of  these  openness  ratios 
is  based  upon  the  assumption  of  100?  covalent  binding.  This  computation 
is  s  impler  than  those  previously  made^  since  the  radius  of  each  atom  doe3 
not  depend  upon  the  other  atoms  present. 

The  information  in  Appendix  A  is  printed  directly  by  the  IBM  computer. 
This  results  in  chemical  formulae  written  in  an  unconventional  form.  Lower 
case  letters  were  not  available.  Hence,  for  example,  the  abbreviation  far 
aluminum  is  written  AL  rather  than  the  customary  Al.  This  will  not  result 
in  confusion  if  one  notes  that  the  abbreviation  for  each  single  element  is 
always  enclosed  between  two  numbers.  As  an  illustration  in  line  1,  one 
finds  the  formula  1AI1B2.  Inasmuch  as  there  is  no  number  between  A  and  L, 
these  letters  designate  a  single  element,  namely  aluminum,  normally  written 
Al.  The  number  1  which  heads  this  formula  indicates  there  is  one  AIB^  group 
in  the  formula.  The  numbers  after  each  element  have  the  meaning  of 
conventional  subscripts.  In  some  cases  rather  unusual  groups,  such  as 
1N1N1  in  line  36,  appear.  The  formula  for  the  compound  is  LiGaN^.  This 
is  written  1N1N1  as  a  matter  of  convenience  in  programming.  The  number 
of  significant  figures  should  be  two  numbers  in  most  cases.  A  larger  number 
are  given  mainly  for  convenience. 


VOLUME  OF  VOLUME  OF  OPENNESS 

CHEMICAL  FORMULA  SOURCE  IONS  FORMULA  WEIGHT  RATIO 

TALI  B2.  *  •  •  7.5073  15.3605  .5113 
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APPENDIX  B 

THE  EFFECT  OF  POROSITY  ON  THE  ELASTIC  PROPERTIES  OF  CERAMICS 

The  sintering  of  pure  ceramic  specimens  generally,  results  in  densities 
lower  than  theoretical  unless  the  particle  size,  forming  pressure,  sintering 
temperature  cycle,  sintering  environment,  and  other  factors  affecting  final 
density  are  carefully  controlled  and  optimized  for  each  new  ceramic  bocfy. 
Lowering  of  density  due  to  the  presence  of  pores  also  lowers  the  elastic 
moduli  values.  In  view  of  the  difficulty  of  obtaining  suitable  ceramic 
specimens  having  maximum  density,  it  appeared  desirable  to  investigate 
methods  of  predicting  the  elastic  properties  of  bodies  having  theoretical 
density  from  measurements  made  on  specimens  with  finite  porosity. 

Mackenzie  J  derived  expressions  for  the  bulk  modulus  (  K  )  and  the 
shear  modulus  (G  )  for  solids  containing  spherical  isolated  pores: 

4K0G0(f~P) 

*  AQo+ZKaP 

G=Gc(f~ZP) 

_  _rfaKo+4Go  \ 

rfxere  ^  ~  ^  +QGoj  (Table  B-l) 

P  ■  volume  fraction  of  pores  -  1  -  .../L- 


Table  D-l 

DIMENSIONLESS  ELASTI C  CONSTANT  OF  OXIDE  CERAMICS 

(USED  IN  POROUS  BODY  COMPUTATIONS) 


CERAMIC  PHASE 

ISJ 

II 

«n 

’3  K0+  1  0o 

9  Ko+  8  3o 

CORUNOUM 

1.97 

MAGNESIUM  OXIDE 

2.03 

SPINEL 

2.01 

FUSED  SILICA 

2.02 

«■-  QUARTZ 

1.98 

PI-1273-M-12 
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Table  B-2 

SHEAR  MODULUS  DEPENDENCE  ON  POROSITY 

BASED  ON  MACKENZIE'S  RELATION 


POROSITY,  P 

Z  -  1.93 

Z  -  1.98 

Z  =*  2.0S 

0 

1.000 

1.000 

1.000 

6* 

0.901 

0.901 

0.898 

10 

0.807 

0.802 

0.797 

20 

0.811 

0.601 

0.681 

»* 

0.123 

0.108 

0.381 

Mackenzie's  model  consisted  of  a  spherical  pore  of  radius  (  Aq)  surrounded 
by  a  spherical  shell  of  solid  with  radius  (  t~a  )  which  in  turn  was  surrounded 
by  an  equivalent  homo  gen  ioua  material  having  the  same  properties  a3  the  pore 
and  shell  combined.  His  solutions  were  obtained  by  expanding  the  function 
ftr  in  a  power  series.  However,  the  higher  power  terms  (i.e.  those 
above  the  third  power)  were  neglected  in  hi3  solution.  Coble  and  Kingery^ 
suggested  the  next  term  in  Mackenzie's  shear  modulus  expression  will  be  of 
the  form  AP2  where  the  constant  ( A  )  may  be  evaluated  by  setting  G  ”  0 
at  P  *  1,  thus* 


G  -Go  (l-ZP+AF?) 


A  -  Z-7 


and  G  =  Go\t -ZP 


Table  B-I  lists  the  values  of  Z  computed  for  several  different  ceramics. 
Since  Z  is  approximately  two  for  many  ceramics,  it  can  be  seen  that  this 
term  increases  the  predicted  value  of  shear  modulus,  slightly,  for  porous 
materials  over  that  obtained  from  Mackenzie* s  original  expression  as  shown 
in  Figure  B-l. 

Coble  and  Kingery^  found  reasonably  good  agreement  between  predictions 
using  the  latter  expression  and  measurements  made  on  specially  prepared 
low-density  alumina  bodies  as  shown  in  Figure  B-l.  However,  recent  data  by 

t  *7 

Lang4*  does  not  fit  either  expression.  In  order  to  investigate  further  the 
effects  of  Mackenzie's  approximate  solution,  an  additional  term  in  the  form 
was  added  to  give 

G  =  G0  (l ~  ZP  -  A  P 2  BP' Z) 

let  G  a  O  at  P  *  1 

and  -”7  *  O  at  P  =  1 

cfP 

then  A  =  2Z  -<3 

B=Z-Z 

G=G0\l-ZPH2Z-3)P2+(2-Z)PS\ 

The  effect  of  the  BPJ  term  was  small  and  the  fit  to  Lang's  data  was 
not  improved.  Clearly,  the  finite  porosity  in  Lang  *s  specimens  was  not  of 
the  dis conti nuous-pore-phase  considered  by  Mackenzie's  model.  In  fact, 
with  the  exception  of  ceramic  bodies  purposely  made  low  in  density  by  foaming 
or  "burn-out  binder"  techniques,  normal  pressing  and  sintering  will  not 
produce  isolated  spherical  pores  except  in  the  final  stages  of  sintering 
where  the  porosity  will  be  under  Mackenzie's  relation  involved  only  the 
elastic  constants  of  the  solid  material  and  the  bulk  density  of  the  porous 
solid.  Hence  it  could  not  distinguish  between  continuous  sol id-phase, and 
continuous  pore-phase  materials. 
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Continuous  pore-phase  materials  could  be  represented  by  a  simple  cubic 
structure  of  equal-size  spherical  grains  in  ■which  an  uncompressed  model  would 
have  a  porosity  of  U8£.  Using  G  =  0  at  P  =  .1*8  for  the  boundary  conditions 
on  Cobles*  extension  of  Mackenzie's  relation  gives 

A~J-.35-Z.09Z 

and  G^Go \~ZP -(4.  35 - 2. 09 Z) P 2] 


It  can  be  seen  that  the  predicted  shear  modulus  is  lowered  by  these  boundary 
conditions  and  approaches  Lang's  data  for  the  continuous  pore-phase  boundary 
conditions.  Therefore,  if  the  constant  A  were  treated  as  an  empirical 
constant  for  continuous  por e-phase  bodies,  Mackenzie's  relation  could  be 
made  to  fit  the  experimental  data.  In  this  case,  the  value  of  A  would 
be  dependent  upon  the  porosity  of  the  unsintered  body  which  is  in  turn 
dependent  upon  the  particle  size,  the  particle  size  distribution,  the 
particle  shape,  the  forming  pressure,  the  percent  binder  and  perhaps  other 
experimental  factors.  Another  approach  to  the  porous  body  problem  was  made 
by  extension  of  relations  derived  by  Kerner  for  elastic  moduli  of  composite 
bodies.  Kerner  derived  an  expression  for  the  shear  modulus  of  composite 
materials  ( Ga  )  contained  within  a  "suspending  fluid"  (subscript  1)  as 


By  considering  a  porous  body  as  composed  of  a  porous  phase  (subscript  2) 
and  a  solid  phase  (subscripts  i  4?  2),  Kerner's  equation  may  be  used  to 
derive  an  expression  for  porous  solids.  For  a  porous  single-phase  solid 
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(  Gz  =  0)  and  the  shear  modulus  becomes: 


Gp  -Cf 


Vi 

15(1 -Vi) 

Q  Vz  Vi 

(7~5iri)G1  15(1 -V,) 


For  comparison  with  Mackenzie’s  expression,  this  equation  may  be  written 
in  terms  of  the  bulk  modulus  (  Ka  )  and  shear  modulus  (  G0  )  of  the  fully- 
dense  body.  Collecting  terms  and  substitution  of 


Z 


f3K0  -4G0 
9Ka  -h  6G0 


gives 


G  -  G, 


1-  P 

’  PZ+(1-P)  ' 


Figure  B-l  illustrates  the  dependence  of  shear  modulus  upon  porosity 
in  alumina  ceramics  according  to  the  two  approaches.  The  later  expression 
is  labeled  "Modify :d  Kerner's  Method"  although  he  did  not  imply  applicability 
of  his  method  to  the  porous  body  case. 

Following  a  similar  line  of  reasoning,  die  bulk  modulus  for  porous 
single-phase  materials  as  derived  from  Kerner's  expressionbecom.es: 

4GqKq  (f-P) 

4G0  +3K0P 

Young's  modulus  for  a  porous  single-phase  materials  was  derived  from 
the  expressions  for  bulk  modulus  and  shear  modulus  and  their  isotropic 
interrelationship: 


5K+G 
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Young’s  modulus  for  porous  media  based  upon  Kerner's  expressions  becomes: 

Co(f-P) 

GO'.ZP _ £. 

Yo  <5  4 

Based  upon  Mackenzie’s  expressions.  Young’s  modulus  becomes: 


Gp 

y0 


GoQ-P)0-ZP) _ 

ZP P  )  / ZPGo  ,  ZP*  ) 

3  4/  [Vo  T  12  j 


A  comparison  of  these  two  expressions  for  the  case  of  porous  alumina  is 
shown  in  Figure  B-2.  The  deviation  of  Mackenzie’s  relation  from  the 
empirical  data  of  Coble  and  Kingery'^  above  30  percent  porosity  is  probably  . 
due  the  higher  order  terms  which  Mackenzie  neglected  diring  his  series 
expansion.  As  Mackenzie  noted,  his  relations  are  valid  only  far  small  porosity. 

For  the  application  of  Kerner's  equations,  the  pore  phase  was  considered 
a3  a  simple  mixture  within  the  solid  phase  without  regard  to  pore  geometry. 

It  might  have  been  expected  that  such  an  approach  would  not  correlate  well 
with  empirical  data  indicated  in  Figures  B-l  and  B-2. 

The  successful  prediction  of  elastic  properties  of  fully-dense  ceramic 
bodies  from  single  crystal  data  described  in  Appendix  C,  postponed  attempts 
to  improve  the  porous -body  model  by  including  terms  dependent  upon  pore 
geometry.  Present  equations  are  reasonably  valid  for  porosities  of  five 
percent  or  less  porosity. 

Young’s  modulus  and  shear  modulus  for  porous  magnesia,  spinel,  mullite, 
and  thoria  are  illustrated  in  Figure  B-3  through  B-8.  For  the  examples  of 
alumina,  magnesia  and  spinel,  values  of  Young's  modulus  >V  and  and 

shear  modulus  Gy  and  GR  for  the  fully-dense  ceramic  body  as  derived  from 
single  crystal  elastic  constants  are  plotted  along  the  ordinate.  In  each 
instance,  the  extrapolated  porous-ceramic  value  lies  between  the  two  values 
computed  from  single  crystal  constants.  This  demonstrates  a  good  correlation 
between  two  independent  methods  for  predicting  elastic  properties  of  fully- 
dense  single-phase  ceramic  bodies.  The  use  of  single  crystal  elastic 
constants  is  described  in  Appendix  C. 
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Figure  &-3  YOUNG'S  MODULUS  OF  POROUS  MAGNESIA 
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Figure  8-4  SHEAR  MODULUS  OF  POROUS  MAGNESIA 
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Figure  B-6  SHEAS  MODULUS  OF  POROUS  SPINEL  (MgO  .  A1 2°3) 
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Figure  B-8  ELASTIC  MODULI  OF  POROUS  THOR i A  (Th02) 


Kerner's  formula  may  also  be  extended  to  the  case  of  porous  multiple- 
phase  ceramic  bodies.  Fcr  example,  the  bulk  modulus  far  a  porous  two -phase 
ceramic  ■would  be  given  by: 


K,  Pi  .  K2P2 
k  +4-Gf)  ysj 3K2 

~  Bl  A  V3 

^G/)  4-G( 

where  \Zj  *  volume  fraction  of  pares  ■  porosity 

subscript  1  refers  to  the  basic  matrix  into  which 
phase  2  material  is  added. 

Application  of  this  expression  to  the  aluminum-silica  system  was  illustrated 
at  the  top  of  Figure  2li.  In  this  system,  an  increase  in  porosity  for  any 
composition  results  in  a  decrease  of  the  composite  bulk  modulus  as  might  be 
expected.  However,  the  mixture  of  phases  with  different  compressibilities 
necessarily  introduces  shear  stresses  between  the  particles  durirg  compres¬ 
sion  of  the  composite  if  the 'mixture  has  been  sintered  or  otherwise  bonded 
together.  In  some  systems  the  presence  of  finite  porosity  will  not  cause 
a  monotonic  decrease  in  the  composite  bulk  modulus  because  of  the  independent 
effects  of  porosity  upon  the  large-scale  shear  modulus.  The  I^yrex-spinel 
system,  for  example,  has  a  compositional  range  in  which  small  porosity 
produces  a  sli^it  increase  in  the  bulk  modulus  (Figure  B-9).  Hiis  rather 
surprising  result  serves  to  illustrate  the  importance  of  including  shear 
stresses  in  any  model  designed  to  describe  the  thermoelastic  properties  of 
composite  media.  It  also  provides  further  evidence  that  Kerner's  method 
for  predicting  thermal  expansion  of  composite  media  should  be  superior  to 
Turner’s  method  in  which  shear  stresses  were  neglected. 

Alternatively,  the  bulk  modulus  of  a  porous  composite  body  might  have 
been  computed  in  two  steps.  First,  the  bulk  modulus  of  the  solid  composite 
could  be  computed  from  Kerner’s  equation.  Then  Mackenzie’s  relation  could 
be  applied  to  the  composite  bulk  modulus  to  obtain  the  values  for  finite 
porosity.  However,  Mackenzie's  relation  always  predicts  a  lowering  of  the 
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moduli  due  to  the  presence  of  pores.  Internal  stresses  arising  from  the 
dissimilar  thermoelastic  properties  of  the  several  phases  would  have  to 
remain  constant  with  increasing  porosity  if  this  approach  were  valid.  But 
in  the  limit  of  high  porosity,  the  internal  stresses  approach  zero.  Hence, 
Mackenzie’s  relation  cannot  be  correctly  applied  to  multiple-phase  ceramics 
or  composite  bodies. 


APPENDIX  C 

RELATIONSHIPS  BETWEEN  SINGLE  CRYSTAL  AND  POLTCRYSTALLINE 
ELASTIC  CONSTANTS  INCLUDING  DATA  FOR  SEVERAL  CERAMIC  PHASES 

The  elastic  properties  of  fully-dense  polycrystalline  bodies  vhich  are 
macroscopically  isotropic  can  be  predicted  from  the  elastic  constants  for 
single  crystals  of  each  of  the  phases  present  by  the  use  of  appropriate 
space  averaging  techniques.  Space  averaging  presumes  random  orientation 
of  the  grains.  The  validity  of  the  random  grain  orientation  equations  will 
be  greatest  if  neither  the  particle  shape  nor  the  forming  technique  induce 
significant  amounts  of  preferred  orientation.  Voigt  took  the  stiffnesses  of 
the  aggregate  to  be  the  space  average  of  the  compliances  of  the  single 
crystal.  Alternatively,  Reuss  took  the  compliances  of  the  aggregate  to  be 
the  space  average  of  the  compliances  of  the  single  crystal.  These  two 
approaches  are  equivalent  to  assuming  uniform  strain  and  uniform  stress 
respectively,  for  each  spherical  grain  in  the  polyciystalline  body,  labile 
neither  condition  will  exist  in  a  real  body,  these  approaches  are  valuable 
in  that  they  define  the  two  extreme  cases.  Empirical  data  for  polycrystalline 
metallic  elements  generally  lie  between  these  extremes  according  to  data 
tabulated  by  Hearmon^f  Similar  results  have  been  found  for  single  -phase 
polycrystalline  ceramics  during  the  course  of  this  study. 

Hooked  law  states  that  stress  (  o'  )  is  proportional  to  strain  (  5  ) 

m 

for  sufficiently  small  strains.  The  generalized  statement  for  an  anisotropic 
medium  may  be  taken  as 

6 

j*f 

vhere  Cjj  *  elastic  stiffness  constants  or  the  moduli  of  elasticity. 
Alternatively,  the  relation  may  be  written  in  terms  of  the  elastic  compliance 
constants  ('  Sy  )  where 

6 

Si  -  Z,  ®ij  *3/ 

.  r1  . 


G-l 
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In  general. 


y 


C/y' 

The  elastic  stiffness  and  complaince  constants  obtained  with  these  relations 

hold  for  a  rectangular  coordinate  system  with  one  preferred  crystallographic 

orientation.  Table  C-I  lists  reported  values  for  oxide  crystals.  If  the 

elastic  constants  for  a  more  general  direction  are  desired,  equations  for 

f>6 

rotated  elastic  constants  must  be  used.  Hearmons  has  tabulated  the  126 
equations  required  for  the  most  general  case.  Ety  integration  of  the  equation 
for  rotated  elastic  constants  over  all  space,  i.e. 


5*y 


2TT 

F(<P,1'iO)d$d{/de 


the  space  averaged  value  will  be  obtained.  This  was  done  by  Hearmon 
who  gives  Voigt's  results  in  terms  of  the  elastic  stiffness  constants  as 

CH  -  j~(3A  +  2B+4C) 

Cr2  = -£-(A+*B-2C) 
q^  =  -L(a-B+ZC) 

where 

•3 A  —  Cfi  * 

•3C  —  C+f  O55+  C„ 

and  the  overline  represents  an  average  value.  Since  space  averaged  values 
are  macroscopically  isotropic,  the  isotropic  constants.  Young's  modulus  (  V  ), 
modulus  of  rigidity  (  G  )  and  bulk  modulus  ( X  )  may  be  derived  in  terms  of 


C-2 
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the  elastic  stiffhess  constants  as  given  by  Rearnon: 

_  (A-8tJC)(A+23) 

2A  +  33+C 


A-B+3C 

5 

A+2B 


u  -  d— r±i_ 

Kv~  v3 

where  the  suffix  v  denotes  the  Voigt  moduli.  Similarly,  Hearaon  derived 
the  egressions  based  on  the  Reuss  approach  as 

KR~  3A+6& 

C  5 

4.A'-4-B'+  3C' 

r  _  5 

**  ~  3Al+2B'+Cr 

where 

3A  -  s,/+  ^22  S33 

Comparison  of  foxing1  s  modulus  or  shear  modulus  as  computed  by  the  Voigt 
and  Reuss  techniques  provides  a  measure  of  the  elastic  anisotropy  of  the 
single  crystal.  For  elastically  isotropic  bodies,  Voigt's  approach  and 
Reuss'  approach  will  give  identical  values. 

For  a  cubic  crystal  such  as  MgO,  only  three  independent  elastic  stiff¬ 
ness  or  compliance  constants  appear 


Cff  -Q>2  -  -‘5-6  *10" dyneejorn2- 

C23  -  Gv  -  C/2  -  87*  IO" dyne,  sjom * 
Co.  =Cs5  -  C/,6  —  14-Q  *10  'dyneojc.fr)2 
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Hence, 


C/t  *•  2Cts 

- 3 - pai 


and  Sit  =  S22-  2^  =  1.08  x  fO'^cm  j/dyne  ' 

-  Sc,/  -  5/2  =-035 xX)'a  crr>*/dyne- 

$44  =s55~  $&= 63b  *10~t9  cuv^ dyne- 


kr=— — 4 —  -  az.z*io6p*i 

*  05/,  +6s,2  r 


In  a  cubic  system. 


Ctf  2.C;2  — 


&„  +2s,2 


and  hence  the  bulk  moduli  (  Ky  and  Kr*  )  will  be  equal  for  all  cubic  crystals. 
Computation  of  both  Kr  and  for  cubic  crystals  provides  a  measure  of  the 
consistency  of  the  reported  values  of  elastic  stiffness  (  Cy  )  and  (  sy  ). 

Birch ^  has  applied  these  space  averaging  techniques  to  several  ceramic 
crystals  to  obtain  the  compressional  wave  and  shear  wave  velocities  for 
polycrystalline  bodies.  For  example,  the  shear  modulus  (  G  )  for  polycrystalline 
corundum  (Al^O^)  computed  from  Birch's  space  averaged  data  arer 

(Voigt)  =•  2l;.0  x  106  psi 

(Reuss)  "  23.2  x  10^  psi 

Langu‘  has  measured  the  dynamic  shear  modulus  of  a  series  of  alumina  ceramic 
bars  with  varying  densities.  Extrapolation  of  his  data  to  theoretical  density 
( yO  =  3*98  g/cnr)  gives  an  experimental  value  of  23. $  x  10  psi  for  the 
shear  modulus  which  lies  between  the  two  predicted  values.  Birch  ^  has 
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tabulated  similar  data  for  eleven  other  mineral  or  synthetic  crystals. 
These  data  have  been  used  to  compute  the  elastic  properties  of  fully  dense 
ceramics  listed  in  Table  XV.  However,  it  is  to  be  noted  that  several  of 
the  ceramic  phases  of  interest  in  studies  of  thermal  expansion,  internal 
stress  due  to  thermal  expansion  anisotropy,  tensile  strength  of  multiphase 
ceramics  and  other  ceramic  properties  have  no  single  crystal  elastic 
constant  data  reported  in  the  literature. 
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